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ABSTRACT 

Results are presented from near-infrared spectroscopic observations of a sample of Bz/T-selected, massive 
star-forming galaxies (sBzKs) at L5 < z < 2.3 that were obtained with OHS/CISCO at the Subaru telescope and 
with SINFONI at the VLT. Among the 28 sBzKs observed. Ha emission was detected in 14 objects, and for 1 1 
of them the [N II]A6583 flux was also measured. Multiwavelength photometry was also used to derive stellar 
masses and extinction parameters, whereas Ha and [N II] emissions have allowed us to estimate star-formation 
rates (SFR), metallicities, ionization mechanisms, and dynamical masses. In order to enforce agreement be- 
tween SFRs from Ha with those derived from rest-frame UV and mid-infrared, additional obscuration for the 
emission lines (that originate in H II regions) was required compared to the extinction derived from the slope of 
the UV continuum. We have also derived the stellar mass-metallicity relation, as well as the relation between 
stellar mass and specific SFR, and compared them to the results in other studies. At a given stellar mass, the 
sBzKs appear to have been already enriched to metallicities close to those of local star-forming galaxies of 
similar mass. The sBzKs presented here tend to have higher metallicities compared to those of UV-selected 
galaxies, indicating that near-infrared selected galaxies tend to be a chemically more evolved population. The 
sBzKs show specific SFRs that are systematically higher, by up to 2 orders of magnitude, compared to those 
of local galaxies of the same mass. The empirical correlations between stellar mass and metallicity, and stellar 
mass and specific SFR are then compared with those of evolutionary population synthesis models constructed 
either with the simple closed-box assumption, or within an infall scenario. Within the assumptions that are 
built-in such models, it appears that a short timescale for the star-formation (~ 100 Myr) and large initial gas 
mass appear to be required if one wants to reproduce both relations simultaneously. 

Subject headings: galaxies: star-formation — galaxies: evolution — galaxies: metallicity — galaxies: high- 
redshift 
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1. INTRODUCTION 

Much action takes place during the ^ 2 billion years of cos- 
mic time between z = 3 and z = 1 .4: both the global star- 
formation rate (SFR) and nuclear (AGN) activity peak at 
z 2, while much of the galaxy assembly and morphologi- 
cal differentiation are taking place and a population of mas- 
sive, passively evolving galaxies is gradually emerging. Thus 
a full exploration of this redshift interval is critical for our un- 
derstanding of galaxy evolution. Unfortunately, this is also a 
relatively difficult redshift range to penetrate observationally, 
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as strong spectral features such as Ca II H&K and the 4000 
A break (for passive galaxies) and [O II] A3727, [O III] A5007, 
Ha and H/3 (for star forming ones) have all moved to the near- 
infrared (NIR). 

Still, most spectroscopic surveys of galaxies at 1 .4 < z < 3 
have used optical spectrographs, relying on much weaker 
spectroscopic features for redshift determination and galaxy 
characterization, with these limitations being partly compen- 
sated by the possib i hty to achieve a fairly high multiplex (e.g. , 
ISteideletal.''2004'; 'Mianoli et al."2005'; 'Le Fevre et al."2005 ; 
Vanzeflaetal. 2005; Lilly et al. 2007; Cimatti et al. 2008; 
Popesso et al. 2009). For passive galaxies at z > 1.4 such 
features practically restrict to a set of absorption features at 
AA ^ 2600-2850 A due to neutral and singly ionized Mg and 
Fe, and for star-forming galaxies to several weak absorption 
lines over the rest-frame UV continuum, most of which due 
to the interstellar medium (ISM) of these galaxies. 

The intrinsic weakness of the absorptions and/or of the 
continuum made such spectroscopic observations very de- 
manding in terms of telescope time. Therefore, several 
studies of large samples of galaxies at 1-4 < z < 2.5 have 
relied on color selections and photometric redshifts. Par- 
tic ularly effecti ve has proven the BzK criterion introduced 
by iDaddi etaH (|2004^, which is able to select both star- 
forming (called sBzKs) as well as passively evolving galaxies 
(pBzKs) over this redshift range. This has enabled estimates 
of SFRs, stellar masses, and clustering properties of such 
galaxies, with samples from ~ lOO's to over ^ 30,000 ob- 
jects (e. g.. iKong et"ani2006l; iDaddi et al.ll2007al; iDunne etal] 
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120091: iMcCracken et alj|2010l) . The BzK technique ensures 
a nearly unbiased selection of z ^ 2 galaxies, including UV- 
selected galaxies and single col or, NIR-sel ected galaxies (e.g., 
iReddy e t al. 2005; McCracken eiiD|2010|). 

Whereas many aspects concerning the evolution of galax- 
ies can be investigated using only photometric redshifts, spec- 
troscopy remains indispensable for a variety of investigations. 
These include full mapping of the local environment (locat- 
ing clusters, groups, filaments and voids), refining SFR and 
mass estimates, measure stellar and ISM metallicities, and fi- 
nally map the internal dynamical workings of galaxies via 3D 
spectroscopy. In this respect, the required telescope time is 
not the only drawback of the optical spectroscopy of galax- 
ies at 1-4 < z < 3. Indeed, optical spectroscopy down to a 
limit as faint as B ^ 25 does not recover but a minor frac- 
tion of the global SFR and stellar mass at z ~ 2, in particular 
missing galaxies that are among the most massive and most 
intensively star-forming ones (Renzini & Daddi 2009). More- 
over, high spatial resolution, 3D spectroscopy of high redshift 
galaxies requires the knowledge of the spectroscopic redshift, 
to make sure that interesting emission lines (e.g.. Ha) are 
free from OH and other atmospheric contaminations (see e.g., 
iGenzel et al. 2006; Forster Schreiber et al. 2009). 

In the case of star-forming isBzK) galaxies, the poor corre- 
lation of mass and SFR with B magnitude is a result of high 
extinction. Therefore, the situation should appear more fa- 
vorable in the NIR, and not only because moving to longer 
wavelengths should reduce the impact of extinction, but also 
because the most active star-forming and most massive galax- 
ies are also among the brightest objects at these longer wave- 
lengths, and one can access strong emission lines such as 
[O II] and Ha. Yet, NIR spectroscopy of z > 1.4 galax- 
ies is still in it s infancy, especially for NIR selected sam- 
ples. lErb et al.] (I2006albircl^ have presented results for over 
100 UV-selected galaxies at z 2, deriving SFRs from the 
strength of Ha. NIR spectroscopic observations of samples 
of z ^ 2 galaxies selec te d in the NIR h ave been also presented 
bv lKriek et all ( l2006allbll2007ll2008allH) . focusing mainly (but 
not exclusively) on passive galaxies by detecting the 4000 
A break, and deriving spectro-photometric redshifts from it. 
NIR spectrosc opy of sBzK galaxies has been carried out by 
iHavashi et alj (l2009i) for a sample of 40 sBzKs, and detected 
Ha emission from 15 of them. Their detections, however, are 
limited to z < 2. Finally, integral field NIR spectroscopy for 
some 60 star forming galaxie s at z 2 has been obtained by 
iForster Schreiber et alJ (|20 09'), for partly UV-selected, partly 
Bz^T-selected targets. Therefore, there is still just scanty spec- 
troscopic information in the rest-frame optical wavelength for 
actively star-forming and heavily obscured galaxies at z — 2, 
many of which would be missed by the UV-selection, and or 
are virtually unreachable by current optical spectroscopy. 

In the perspective of improving upon this situation, in 2004 
we started NIR spectroscopic observations of z ^ 2 galax- 
ies primarily selected on BzK technique, and using a vari- 
ety of NIR instruments, namely OHS, CISCO and MOIRCS, 
at the Subaru telescope and SINFONI at the VLT. Our in- 
tent was to explore the effectiveness of NIR spectroscopy to 
improve our characterization of z ^ 2 galaxies using a rela- 
tively small pilot sample of them, while assessing the feasi- 
bility of wider surveys with future instruments with higher 
multiplex. In this paper we present the results of observa- 
tions with the OHS/CISCO instruments on Subaru and SIN- 
FONI instrument on VLT of sBzK galaxies, leaving the re- 



sults obtained with the MOIRCS instrument for a future pa- 
per With the observations presented here we have attempted 
to measure for each galaxy the dust extinction, SFR, ISM 
metallicity, and dynamical mass, while checking for a pos- 
sible AGN contribution. Physical quantities are derived as- 
suming the concordance cosmology, i.e., JIm = 0.3, fl\ = 0.7, 
and Ho = 70 km s~' Mpc~^ and photometric magnitudes are 
expressed in the AB system jOke & Gunnll983i) if it is not ex- 
plicitly noted otherwise. F or the solar oxygen abunda nce, we 
use 12 + log(O/H)0 = 8.69 (Allende Prieto et al. 2001). Emis- 
sion line width are measured assuming a Gaussian profile 
and the FWHM and the velocity dispersion (a) is given by 
FWHM = 2.355(7. 

2. SAMPLE SELECTION, OBSERVATIONS, AND DATA 
REDUCTION 

2.1. Sample Selection 

The sBzK galaxies which are the object of the present study 
have been culled from the ^T-band selected catalog of objects 
in the EIS Deep3a field and Daddi field, which a re des cribed 
alon g with the photometric data in iKong et aTl (120061 here- 
after lKOd) . We have primarily selected sBzKs with spectro- 
scopically confirmed redshifts from the rest-frame UV spec- 
troscopy obtained with the VIMOS instrument at the VLT 
(Daddi et al., in preparation). Then we have selected those 
for which the Ha emission line falls at wavelengths where 
atmospheric and instrumental transmission are high and is 
not contaminated by OH-airglow lines. Additional criteria 
have then been applied to select objects which (1) are bright 
in the /T-band (Kab ^ 22) to select massive galaxies, i.e., 
> 5 X 10"' Mq (iDaddi et al.l l2004ai see also EB, (2) 
have red z-K colors, (z-K)ab ^ 2 to select possibly more 
reddened ones, (3) are not detected in X-rays down to (2- 
7) X 10~'^ ergs s"' cm"^ (if in the Daddi field) and do not 
show strong AGN features in rest-frame UV spectra, so to ex- 
clude galaxies with an obvious AGN contribution, and (4) are 
bright in Spitzer/MlPS 24 pm images, with flux > 70 /iJy (for 
the Daddi field) and > 100 /iJy (for the Deep3a field), to make 
sure that they are actively forming stars. 

Our selected objects in the Daddi field are not detected in 
XMM-NewtonI Chandra X-ray data with 0.5-2 keV luminosi- 
ties Lx ^ lO'*^ ergs s"', which ensures the exclusion of lumi- 
nous unobscured AGNs (M. Brusa, private communications). 
However, one sBzK in our sample shows mid-IR (MIR) ex- 
cess which suggests it may host an obscured AGN (see ^5.31 
Daddi et al. 2007a b). We note that only 10 of our 28 objects 
comply to all the above constraints, and pairs of sBzKs that 
can be put in the same slit have been preferentially included 
to have a larger sample. Among the 28 objects selected for 
NIR spectroscopy, 20 objects have a redshift from rest-frame 
UV spectra, and 13 objects satisfy all 4 criteria above and 
have a spectroscopic redshift. Object IDs, coordinates, red- 
shifts from rest-frame UV spectra, and magnitudes are listed 
in TablelU 

2.2. OHS/CISCO Spectra 

Among the selected objects, 7 sBzKs were observed on May 
1, 2004 and April 24-26, 2005, by using the OH-airglow sup- 
pressor (OHS; Iwamuro et al. 2001) with the Cooled Infrared 
Came ra and Spectrograph for OHS (CISCO; iMotohara et alj 
120021) mounted on the Nasmyth focus on the Subaru tele- 
scope. The JH grism covering the 1 . 1-1 .8 pm spectral range 
and a 1" slit were used, which give a spectral resolution 
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R ~ 200. On May 6, 2004 and April 30, 2005 K-hmd spec- 
troscopy with the wK grism covering the 1 .85-2.5/im spectral 
range was carried out for 5 sBzKs with CISCO (i.e., with- 
out OHS), which with a slit width of 1" gives a resolution 
R 2± 300. Position angles were selected to put bright nearby 
objects into the slit, so to facilitate target acquisition since our 
targets are usually too faint to be seen on acquisition images. 
In case no nearby object was available, target acquisition was 
made with blind offset from bright objects as close as pos- 
sible to the target object. Each spectrum was obtained by 
taking 4 x (900s or 1000s) exposure sequences with slightly 
modified ABBA standard nodding pattern to avoid bad pixels. 
As spectroscopic standards to correct for atmospheric and in- 
strumental transmission, the A-type stars SAO 120721, SAO 
121153, SAO 122123, SAO 180911, and SAO 180521 were 
observed for the JH grism, while the white dwarf GD 153 and 
the A-type star SAO 121856 were observed for the wK grism. 

All the data were reduced with standard procedures. Two 
dimensional spectra were produced with flat fielding, distor- 
tion correction, bad pixel rejection, and sky subtraction by 
making a median sky from adjacent detector areas. Sky resid- 
uals were subtracted by fitting polynomials in the spatial di- 
rection, and then co-addition was carried out with appropriate 
offsets for the dithering width. The resulting 2D spectra were 
then collapsed to ID spectra. The tilt of the spectra on the 
array was corrected by adopting the tilt of standard star spec- 
tra. The wavelength calibrations were carried out by using the 
standard pixel-wavelength relation of OHS/CISCO with sys- 
tematic error of < 0.5 pixels (3A an d 5A for OHS and CISCO, 
respectivelv: lMotohara et ani2001h . Atmospheric and instru- 
mental transmissions were corrected by using the ID spectra 
of the standard stars reduced in the same way as object frames. 

The absolute flux calibrations were carried out by compar- 
ing the photometric J- and /T-band fluxes with those derived 
from the object spectra convolved with the filter transmission 
curves. Here we have adopted 2" aperture magnitudes, hence 
assuming that in the observed spectra both continuum and 
lines come from same region of a galaxy. Noise spectra were 
derived by measuring rms of counts in the blank sky region. 

2.3. SINFONI Spectra 

On April 14-16, 2005 we observed 16 sBzKs with the Spec- 
trograph for Integral Field O bservations in the Near Infrared 
(SINFONI; EisenhaueretaD|2003; Bonnet et al. 2004) at the 
VLT/UT4. The grating covering the 1. 45-2 .45 /im spec- 
tral range with R ~ 1500 and a pre-optics giving a 125 x 250 
mas pix~' spatial resolution were used. Each object was ob- 
served by 1 or 2 sequence(s) of 4 x (450s or 900s) exposures 
with 4 arcsec dithering both in x- and y-direction. The B-type 
stars, HIP 007873, HIP 026816, HIP 031768, HIP 068100, 
HIP 068372, HIP 072367, HIP 075711, HIP 083861, HIP 
092957, HIP 094333, HIP 095806, and HIP 099244, were 
observed as spectroscopic standards. 

The reductions of SINFONI data were carried out with 
the SINFONI pipeline" and custom scripts, including flat- 
fielding, sky-subtraction by median sky of a sequence, bad 
pixel rejection, distortion correction, wavelength calibration 
by arc lamp frames, cube reconstruction, residual sky subtrac- 
tion by subtracting the median in x-direction from each pixels, 
and finally co-addition with appropriate offsets. Telluric and 
instrumental absorptions were corrected with spectroscopic 

' ' http://www.eso.org/sci/data-processing/software/pipelines/ 



Standard star frames. Since neither the standard stars nor the 
objects are affected by slit losses, the absolute flux calibra- 
tions were obtained by dividing the object spectrum by the 
standard star spectrum, scaled to the broad band photometric 
magnitudes. 

3. EMISSION LINE MEASUREMENTS 

Among the 28 sBzKs observed, emission lines were de- 
tected in 13 of them, and identified as Ha in all cases. Emis- 
sion lines in the SINFONI data cube have also been detected 
in an object in the vicinity of Dad-2426 (hence named Dad- 
2426b), and the lines are identified as Ha and [N II]A6583 at 
z = 1 .772 according to the separation between them. Optical 
composite images and /T-band images of the galaxies with Ha 
detections are shown in Figure [T] and their spectra around Ha 
are shown in Figures|2]and[3]for OHS/CISCO and SINFONI, 
respectively. 

Due to their low resolution in wavelength, the Ha 
emission is always blended with [N II]AA6548,6583 in 
the OHS/CISCO spectra, while these lines are resolved 
in the SINFONI spectra. We have measured the Ha 
line fluxes by fitting multiple-Gaussians for Ha and 
[N II]AA6548,6583, assuming that all three lines have same 
width, [N II]A6583/[N II]A6548 = 3, and a constant contin- 
uum flux. This leaves redshift. Ha emission line flux /(Ha), 
flux ratio /(Ha)//([N II]A6583), line width a, and constant 
continuum flux as free parameters. For most of the objects, 
the fitting converged to a single solution against various initial 
guesses, and the derived line widths agree with those expected 
from instrumental resolution in the case of OHS/CISCO spec- 
tra. However, there are several exceptions. For D3a-8608 and 
Dad-2426 the procedure did not converge to stable solutions, 
and we fixed the line width to that of the instrumental pro- 
file (measured on the OH lines close to the position of emis- 
sion lines). Small perturbations to the assumed line widths 
did not change fluxes and equivalent widths (EWs) apprecia- 
bly. Since the [N II]A6583 line is below the detection limit 
in D3a-3287 and D3a-4626, only a single Gaussian to the Ha 
line was fitted. D3a- 11391 appears to have both broad and 
narrow Ha compon ents. This object also shows MIR-excess 
jDaddi et al ."2007b') as discussed later ( ^5.3t , hence the spec- 
trum was fitted with both broad and narrow Ha lines, which 
resulted in a better compared to the fit with only narrow 
line components. We could not find a good fit when includ- 
ing the [N II] lines, possibly due to residual of OH sky lines 
and rapidly variable atmospheric transmission at the edge of 
//-band. The best fit Gaussian functions are overplotted in 
Figure 12] and [3] 

The resulting redshifts derived from Ha, along with the 
fluxes and EWs of Ha and [N II]A6583, and velocity disper- 
sions are listed in Table [3] The EWs listed in Table |3] are 
derived with the measured line fluxes and continuum fluxes 
from the best-fit spectral energy distribution (SED) (see ^4.1t . 
Ha EWs are also corrected for stellar Ha absorption, derived 
from the synthetic spectrum that best-fits the SED. The error 
in the continuum flux from the best-fit SED around the posi- 
tion of Ha is estimated to be 20%. For Dad-2426b, which 
is not listed in our /T-selected catalog, the continuum flux 
is calculated from its 5<t limit of K^b =21.5, and the stel- 
lar Ha EW is assumed to be 4.4 A, which is the average 
of stellar Ha EWs estimated for the other objects. These 
redshifts agree well with those from rest-frame UV spectra, 
except for 2 objects, Dad-2426 (zuv = 2.36 and zhu = 2.40) 
and D3a-6397 (zuv = 2.00 and zho = 1.51). The spectrum of 
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Dad-2426 is located a little far from the detector center and 
since detector distortion increases with distance from the cen- 
ter, we first doubted that the adopted pixel-wavelength rela- 
tion may not apply in such case. However, a cross-check of 
wavelength with OH lines showed that the wavelength was 
well calibrated. This object was also observed with SINFONI 
and a feature is marginally seen at 2.23//m which if due to 
Ha corresponds to z = 2.40, consistent with the redshift from 
the CISCO spectrum. For D3a-6397, the discrepancy is very 
large, but its zuv is quite uncertain whereas [N II]A6583 is 
well detected at the expected position relative to Ha, hence 
we consider zho as more reliable. 

From the Hq;/[N II]A6583 emission line ratio it is pos- 
sible to estimate whether there is a significant contribu- 
tion from an AGN component, though more emission lines 
such as H/3 an d [O ml are required for a more robust di- 
agnos t ic (e.g., [V eilleux & Osterbrock' '19871; [Baldwin et alJ 
119811: iKauffmann et al. 2003). Here we classify sBzKs 
with [N II]A6583/Ha < 0.7 as star-formation dominated 
and those with [Nil] A 65 8 3 /Ha > 0.7 as AGN dominated 
jSwinbank et alj|2004l) . Adopting this criterion, one sBzK, 
D3a-8608, is classified as AGN-dominated. Besides the emis- 
sion line ratio diagnostics, D3a-1 1391 could host an AGN be- 
cause of its broad-line component. This would read an AGN 
fraction of ~ 15% for our sample of 14 galaxies with de- 
tected emission lines. Note that true AGN fraction among 
sBzKs might be even higher, since part of the objects have 
been pre-selected for lacking AGN features. In any case, 
this AGN fraction is not far from that (~ 25%) estimated 
for th e full sBzK population in the Deep-3a and Daddi fields 
(lK06h . while up to ~ 50% of massive (M* ~ IO^Mq) sBzKs 
show MIR-excess, possibly related to heavily obscured AGN 



show MlK-excess, p c 
jDaddi et al.ll2007bf) . 



4. PHYSICAL PROPERTIES OF sBzK GALAXIES 
4.1. Stellar Masses and Reddening 

Stellar masses (M*) and the reddening E(B-V) were de- 
rived by SED fitting to the BRIzJK- and BRIK-banA data for 
the sBzKs in the Deep3a field and the Daddi field, respec- 
tively. Spectroscopic redshifts determined from Ha emission 
lines were used for the fit. 

Template SE Ds were generated by using PEGASE.2 
jpioc & Rocca-Volmerange 1 997^ 19^ ), assuming a Salpeter 
initial mass function (IMF; iSalpetei? 119551 for stars with 
O.IMq to 12OM0. These models assume exponentially de- 
clining SFRs and incorporate chemica l evolution using chem- 
ical yields from 'Woosl ev & WeaveJ ([T995). Real galax- 
ies at z ~ 2 are unlikel y to have evo lved with exponen- 
tially dechning SFRs (cf. ICimatti et alj '2008; Renzini 2009; 
iMaraston et alJlToiOl) . and their gas accretion histories may 
be radically different compared to those assumed in the PE- 
GASE.2 models. Although the use of these models intro- 
duces some rigidity in the SED fits, we believe that the de- 
rived stellar masse s should be correct within a factor of ^ 2 
(iDrory et alJl2004l) . A systematic overestimate by up a factor 
~ 3 may be present for those galaxies in which the bulk of 
stars have ages around ^ 1 Gyr, i.e., at the peak of the contri- 
bution by the TP-AGB phase of stellar evolution, that was not 
adequately included in the PEGASE.2 (e.g.,LM araston 1998, 
l2005HKaiisawa et alj2009HMagdis et all2009 ). Furthermore, 
a systematic reduction by a factor ~ 1.6-2 of both SFRs and 
masses would be produced adopting a bottom light IMF, such 
as those of .Kroupa (.2002) or .Chabrier (.2003.) . 



The adopted models assume SFR e-folding times from 
Tsf =100 Myr to 500 Gyr with various age grids ranging from 
10 Myr to 15 Gyr. The synthetic SEDs are then dust at- 
tenuated according to the Calzetti law (Calzetti 2001) with 
< E(B-V) < 1.5. Absorption by neutral hydrogen inter- 
vening along the li ne of sight is also applied (Madau 1995]; 
iMadau et"an Il996h . Fitting is finally obtained by a x^- 
minimization, with fixed redshift derived from Ha, and with 
the constraint that age cannot exceed the age of the universe at 
the observed redshift. Stellar masses are then derived by mul- 
tiplying the luminosity by the M/L ratio of the model which 
best-fits the observed SED. 

The best-fit SEDs together with observed data points are 
shown in Figure |4] and the derived stellar masses and redden- 
ings are listed in Table |4] For a consistency check, the stellar 
masses and reddenings are also calculated by using t he BzK- 
calibrations as given by Equation (4) and (6) in iDaddi et alJ 
(l2004ai) : 

E(B-V) = 0.25(B-z+0.1)ab, (1) 



log(M,/lO"M0) = -0.4(/:i°^-21.38) 

+0.218[(z-/:)ab-2.29], 



(2) 



which are based on SED fits to the full UBVRIzJHK -hand 
data (cf. BRIzJK in this study). These BzK-hased stellar 
masses and reddening parameters are also reported in Table 
[T] Figure |5] shows the difference of stellar mass and red- 
dening between the BzK-hased values and those from SED 
fitting, for the Ha detected objects, as a function of stel- 
lar masses from SED fitting. Average offsets and disper- 
sion between two estimators are {log{M^,5ED / M^^ezk)) =0.03 
and o'(log(M*sED/-W*B a:)) = 0.22, respectively, for the stellar 
masses, and {AE(B-V)) = -0.03 and (7(E(B-V)) = 0.07, re- 
spectively, for the reddening. Hence, they agree reasonably 
well and without large systematic offsets. 

4.2. Ha Star-formation Rates 

The Ha luminosities are converted into SFRs following the 
relation in lKennicut3 (119981) : 



SFR(MQyr-') = 7.9 x IQ-^^L^Jtrgs. s"'), 



(3) 



which assumes solar abundance, a Salpeter IMF, and constant 
SFR within the last 100 Myr. Although the conversion factor 
may vary from 2.6 x 10""*^ to 8.7 x 10""*^, depending on metal- 
licity, IMF, and star-formation history (Buat et al. 2002), we 
adopt Kennicutt's conversion as it is the most commonly used 
and it simplifies the comparison with results from the litera- 
ture. 

Although Ha is usually considered to be relatively un- 
affected by dust extinction, dust extinction must be taken 
into account for galaxies like sBzKs which undergo vigor- 
ous star-forma tion and are hkely to be heavily obscu red by 
dust (Table l4l IDaddi et alJl2005t iPannella et al.ll2009l) . The 
H/3 emission line would enable us to estimate the dust extinc- 
tion via the Balmer decrement, but unfortunately it is not de- 
tected in our spectra. The 3a upper limits for the H/3 flux 
of Ha detected sBzKs is (Ha/H/3) 

limit > 2.86, which can- 
not significantly constrain the reddening, except for D3a-475 1 
whose Ha /H/3 > 5.0 implies stellar E{B-V) > 0.2, consis- 
tent with the stellar E(B-V) = 0.25 from SED fitting. Corre- 
spondingly, the reddening values from SED fitting are used 
for the extinction correction of the Ha luminosities. One 
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may expect that nebular emission lines from H II regions 
could suffer from larger dust e xtinction than the stellar con- 
tmuum enussion ('Calzetti"2001). To cope with this problem, 
|Cid Fernandes et al. (2005) compared the extinction of stel- 
lar continuum with that of nebular emission by using star- 
forming galaxies from the SDSS data set, and found a linear 
correlation between them. Thi s relation was combined with 
the Calzetti extinction curve bv lSavaglio et alj ( l2005i) . obtain- 
ing: 

Ay = 3.173+ 1.841A*-6.4181og(^^^ , (4) 

where Ay and Ay are gas and stellar visual extinctions, re- 
spectively, and (HQ!/H/3)th is the line flux ratio calculated 
from the atomic physics theory. Here we assumed the case 
B recombination with T = lOOOOK and (Ha/H;9),h = 2.86 
jOsterbrock & Ferlandll2005l) . Extinction coiTected Ha SFRs 
are then derived by using Ay from the above relation. Vi- 
sual extinctions of stars and gas, extinction at Ha, Ha lumi- 
nosities, and SFRs with and without extinction corrections are 
listed in Table |5] 

After extinction correction, Dad-2426 and D3a-6397 show 
SFR > lOOOMeyr"'. Dad-2426 has been detected at 1.2 mm 
at more than 3a level, and its SFR from FIR luminosity 
of 8OO-19OOM0yr-' (iDannerbauer et al.ll2006l) is consistent 
with the Ha-based SFR. In our sample with Ha detection, 
more than half of the sBzKs have SFR > lOOMoyr"', after 
extinction correction. These Ha SFRs > lOOMoyr"' agree 
with the SFRs derived by the combination of UV and MIR 
for the sBz Ks w ith similar range of stellar mass in GOODS 
fields (see ^5.31 Daddi et al. 2005, 2007a). The comparison 
between SFRs from different indicators and diff erent extinc- 
tion corrections will be further discussed in ^5.31 

4.3. Metallicities 

Accurate metallicities for the ionized gas can be ob- 
tained once the electron temperature is derived from 
the ratio of auroral to nebular emission lines, such as 
[O III]AA4959,5007/[O III]A4363. However, auroral lines 
are intrinsically faint, in particular in the metal rich galax- 
ies, since the electron temperature decreases due to efficient 
cooling by metal lines. Hence, detection of such lines in 
high-z galaxies is not feasibl e with current facilitie s except 
for strongly lensed galaxies (I Yuan & Kewlevll2009h . Alter- 
natively, abundance indicators using strong emission lines are 
widely used to determine the metallicity, being calibrated with 
the Te method and/or photoionizatio n models. One of the most 
commonly used indicators is /?23 ( iPagel et al.lll979h . which 
is defined as XogRn = log[([0 II]A3727 + [O III]A4959 + 
[O III]A5007)/H/3]. However, /?23 cannot be applied to our 
sample either, due to the lack of required emission lines. 
Here we use instead the N2 index (St orchi-Bergmann et al.l 
11994ft defined as N2 = log([N II]A6583/Ha). Since Ha and 
[N II]A6583 are close to each other, the N2 index has the ad- 
vantage of being insensitive to dust extinction and flux cali- 
bration, though the origin of nitrogen is rather complicated. 

Metallicities derived from different calibrations are known 
to show (almost systematic) discrepancies in the mass- 
metallicity relation of the lo cal s tar-forming galaxies 
dTremontietal.' '200 4 hereafter IT04I) of up to 0.5 dex 
llllison 2006; Kewle v & Ellisonl2008l) . Thus, it is essential to 
derive metallicities in a consistent way, as later we will com- 
pare the metallicities of our sBzK sample with those of galax- 



ies at different redshifts drawn from the hterature. Bas ed on 
iKewlev & Dopltal ( l2002l) . lKobuinickv & KewlevI dlOOl here- 
after KK04) derived an analytical expression for the conver- 
sion of N2 into 12 + log(0/H) that is consistent with the /?23 
calibration. The relation is expressed as: 

12 + log(0/H) = 7.04 + 5.28Xn„+6.28X4+2.37Z4 
-log^(-2.44-2.01XN„ 
-0.325X4 + 0.128X4) 
+ 1 O^Nii-0.2 log ^(-3 . 1 6 + 4 .65Xn„) , 

(5) 

where Xn,, = logEW([N II]A6583)/EW(Ha) and q is the ion- 
ization parameter. Since the EWs of [N II] and Ha were de- 
rived by assuming a flat continuum, EWs in can be re- 
placed by the corresponding emission line fluxes. According 
to KK04, the ionization parameter is derived iteratively from 
the ionization-sensitive index log = log[([0 III]A4959 + 
[OIII]A5007)/[OII]A3727]. However, once more also O32 
cannot be derived for galaxies in our sample, as some of the 
required lines are not available. We have then assumed a 
constant ionization parameter, an assumption that is justified 
as follows after comparing metallicities from the N2 index 
and those from the RnjOyi index. To derive metallicities 
from both indices, we used emission line fluxes for 75,561 
star-forming galaxies in the SDSS DR4 archive compiled by 
MPA/JHUcollaboration'2. Then we calculated 12+log(0/H) 
for these SDSS star-fo rming galaxies using Equation (|5]) and 
by the relation (IKK04 : 

12 + log(0/H)= 9.11-0.218x-0.0587;i:2 
-0.330x^-0.199/ 
-y(0.00235-0.01105jt:-0.051x2 
-0.04085x^-0.003585/), (6) 

where x = log(/?23) and y = log((932). 

Figure |6] shows the correlation between these two metal- 
licity calibrators, with different ionization parameters q = 
1 X 10^ 2 X 10\ 3 X 10^ and 4 x 10''. This figure indicates 
that a value q = 3 x lO'' makes the two metallicities to agree 
for galaxies near the peak of the metallicity distribution, while 
the N2 calibration tends to overestimate the metallicity for 
12 + log(0/H)R,3 < 8.8. Given that the metallicities of our 
sample are generally large, 12 + log(0/H) > 8.8 (see below), 
we have adopted ^ = 3 x 10^ in Equation (|5]l. 

The resulting gas-phase oxygen abundances of our sBzKs 
are Hsted in Table|6] Although D3a-8608 has [N II]/Ha > 0.7, 
which indicates that the ion ization may be dominated by an 
AGN (ISwinbank et alJl2004h . we derived its metallicity with 
the same equation, just for a reference. The average value of 
12 + log(0/H) in our sample is 8.97±0.21 excluding objects 
with upper limits and AGN dominated features, i.e., D3a- 
3287, D3a-4626, D3a-8608, and D3a- 11391. 

4.4. Dynamical Masses 

Thanks to the higher spectral resolution (R ~ 1500) of 
the SINFONI instrument. Ha and [N II] emission lines 
are well resolved, which enables us to measure the veloc- 
ity width of individual objects. In the following we refer 

http://www.mpa-garching.mpg.de/SDSS/ 
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to such velocity width as the velocity dispersion, although 
we are aware that in many cases it is likely due mostly 
to ordered rotation rather th an true velocity dispersion (cf. 
iForster Schreiber et alJl2009l) . The measured velocity disper- 
sions, after deconvolving in quadrature the instrumental res- 
olution (R = 1500, or (Tinstrument = 85 km s"'), are listed in Ta- 
ble |3] TTien_dYnaniical masses a re derived using th e equation 
(e.g., iBinnev & TrernainQil987HPettini et alj|2001h : 



where for the half-Ught radius, r^i we have adopted 6 kpc, i.e., 
the average value for th e sBzKs in the K20 survey measured 
bv lDaddietaT](l2004bl) on the HST/ACS F850LP-band im- 
age. This assumes that the Ha emission comes from H II re- 
gions which are virialized within the potential well of the host 
galaxy, and that the spatial distribution of H II regions repre- 
sents that of underlying light distribution. Dynamical masses 
derived in this way are listed in Table|7] Note that the velocity 
dispersion of D3a-4751 is less than the width of the instru- 
mental profile, hence the velocity dispersion could have large 
uncertainty. The average of dynamical masses of 8 sBzKs 
with measured velocity dispersion is (Mdyn) = 2.3 x 10 "Mq 
with la scatter of 1.5 x lO^M©. 

There are several possible sources of error in these esti- 
mates of the dynamical mass. The half-light radius which 
is assumed as constant must actually vary from one object to 
another. For examp le, in the cas e of the UV-selected galax- 
ies studied by Erb e t al.l (l2006bh the half-light radius varies 
by more than a factor of 5 from, ^ 2 kpc to ^ 11 kpc while 
the average is ~ 6 kpc. In a re cent integral field Ha spectro- 
scopic survey of z ~ 2 galaxies lForster Schreiber et al.l (l2009h 
find an average half-light radius of the Ha emitting region of 
3.4 kpc. If we adopt this value, the derived dynamical masses 
would be 0.6 times smaller than estimated above. The numer- 
ical factor in Equation (|7]l could also be a dominant source 
of uncertainty. The adopted value of 5 is valid for a sphere 
of uniform density, while the actual value of this factor de- 
pends on various parameters, such as the mass distribution 
within the galaxy, its velocity structure, and the relative con- 
tributi ons of rotation and press ure support for star-forming re- 
gions (Lanzoni & Ciotti 2003). Considering a disk geometry, 
|Erb et al. (2006b) used 3.4 for the factor and obtained the av- 
erage dynamical mass of UV-selected star-forming galaxies at 
z~2ofMdyn 1 X 1O"M0 with Ict scatter of 8.5 x IO'^Mq, 
after accounting for differences in scaling factor As men- 
tioned above, several of these sBzKs are likely to be rotating 
disks (e.g., Genzel et al. 2006), in which case the reported dy- 
namical masses are overestimated by ^ 50%. Therefore, the 
dynamical masses derived here could be taken as upper limits 
unless half-light radius is larger than 6 kpc. 

For most of the objects these dynamical masses are ^ 2-3 
times larger than the stellar masses, which in pri nciple could 
be due to th e presence of large gas masses (e.g.. iDaddi et al.l 
120081 120T0I) . However, spatially resolved integral field spec- 
troscopy would be required for a more robust determin ation 
of the dynamical mass (e.g. . IForster Schreiber et al]|2009i) . 

4.5. Composite Spectrum ofsBzKs 

The composite spectrum shown in Figure [T] is made by 
stacking the 6 SINFONI spectra of sBzKs with Ha detec- 
tion and without AGN features, where Ha and [N II]A6583 
are well resolved. Note that we have excluded D3a- 11391 



which shows broad-line Ha component as well as the sBzKs 
with higher [N II]/Ha ratio as AGN candidates. Individual 
spectra were corrected to the rest-frame wavelength and then 
stacked with weights inversely proportional to the square of 
the IfT rms versus wavelength. Motivated by the finding of a 
broad-line Ha component in the stac ked spectra of z ^ 2 star- 
forming galaxies in the SINS survey jShapiro et al.ll2009b . we 
also fit multiple Gaussian functions to the composite spec- 
trum, with and without a broad Ha component, though due to 
the shorter exposure time for each object and smaller number 
of spectra used for the stacking, the S/N is lower than that of 
the SINS stacked spectrum. Including broad a Ha component 
makes the reduced-x^ of the fit about 25% smaller The best 
fit Gaussian profiles are also shown in Figure|7] 

The spectral properties derived from this composite spec- 
trum are EW(Ha) = 102 A, EW([N II]A6583) = 41 A, and 
FWHM = 350 km s"' or a = 150 km s'^ for the fitting 
without a broad Ha component, and EW(Ha; narrow) = 
50 A, EW(Ha; broad) = 54 A, EW([N II]A6583) = 24 A, 
FWHM(narrow) =170 km s"' or (T(narrow) = 73 km s"' 
and FWHM(broad) = 870 km s"' or cr(broad) = 370 km s"' 
for the fitting with a broad Ha component. Velocity dis- 
persions above, are corrected for the instrumental resolution. 
Gas-phase oxygen abundances derived from the narrow line 
components of the composite spectrum can be calculated as 
12 + log(0/H) = 9.12 and 9.03; the line widths lead to a dy- 
namical mass of Mdyn = 3.7 x 1O'"M0 and 1.5 x lO^M© for 
the fittings with and without a broad component, respectively. 
If the broad component is real, the dynamical mass is a fac- 
tor ^ 2 lower than the average stellar mass of the galaxies 
used for the stack (7.6 x lO'^M©). However, stellar masses 
have been derived here assuming a straight Salpeter IMF, with 
most of the mass being provided by low-mass stars. Turn- 
ing from a bottom-heavy IMF, to a bottom-light one such as 
the IMF of Chabrier (2003), the average stellar mass drops to 
4.5 X 1O'"M0, quite close to the estimated dynamical mass, 
and we note that a bottom-light IMF appears to be more ap- 
propriate to accou nt for the obse rved mass to light ratio of 
local galaxies (e.g.. lRenzinill2005l) . In addition, we may have 
somewhat overestimated stellar masses as our SED fits are 
based on stellar population models that do not incorporate 
the TP-AGB phase of stellar evolution. We conclude that the 
discrepancy between our dynamical and stellar masses is pri- 
marily a result of the assumed IMF, and regard the dynami- 
cal mass derived including the broad Ha component as quite 
plausible and close to the actual stellar mass as well. 

5. DISCUSSION 

5.1. sBzKs with and without Ha Detection 

In most cases, non-detection of Ha could well be due to a 
combination of bad weather, poor seeing, mis-alignment dur- 
ing the blind offset, and bad transmission or strong OH emis- 
sion at the wavelength of the lines. However, it is interesting 
to see whether there are systematic differences in colors and 
physical properties between sBzKs with and without Ha de- 
tection. In Figure[8]the (B - z)-(z- K) color-color diagram for 
the objects in our /T-selected catalog is shown, using different 
symbols for Ha detected and non-detected objects. Ha de- 
tected sBzKs have on average ^ 0.15 mag bluer BzK colors 
than those without Ha detection, which indicates a systemat- 
ically lower dust extinction because the reddening vector runs 
parallel to the diagonal line separating sBzKs from the other 
objects (.Daddi et al...2004a) . Ha detected sBzKs also appears 
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to be slightly brighter in K-hund compared to those without 
Ha detections. Therefore, Ha detections seem to be biased in 
favor of more massive, but less extincted objects. 

The reddening E{B - V) from Equation ([T]i and the ex- 
tinction corrected SFRs derived from the flux in the rest- 
frame UV (observed- frame B-band) according to the recipe 
in Da ddi et al.l (l2004ah are plotted in Figure|9]as a function of 
stellar masses derived from Equation (|2]i. There are no sig- 
nificant differences in stellar masses (~ 0.1 dex) and SFRs 
(20%, but uncertainties are large) between Ha detections and 
non-detections, while the average reddening is 0. 1 mag higher 
for the Ha non-detections. Therefore, the primary factor af- 
fecting the detectability of Ha emission lines from sBzKs ap- 
pears to be the amount of dust extinction. This is also seen in 
iHavashi et alj ( 1200% where most of sBzKs with non-detected 
emission lines have redder BzK colors and large E(B-V) val- 
ues, > 0.5. 

5.2. Comparison with Other Z — 2 Ha Spectroscopic Survey 

Figure [TO] compares ^T-band magnitudes, stellar masses, 
reddenings. Ha fluxes and Ha luminosities of our BzK sam- 
ple with those from other Ha spectroscopic surveys at z ~ 2. 
The Ha fluxes and luminosities are not corrected for ex- 
tinction. The galaxy sa mples include the rest- frame UV- 
selected BX/BM galaxies (|Erb et al."2006b'c'), rest-frame op- 
tically selected sBzK galaxies (Hayashi et al. 2009), and the 
z — 2 star-forming galaxies observed by the SINS survey 
jForste r Schreiber et al. 2009) which includes BX/BM galax- 
ies as well as Bz/T-selected galaxies. Compared to the other 
Ha detected star-forming galaxies, our sample contains only 
/^-bright, massive objects mainly because of the brighter K- 
band limiting magnitud e of the imaging survey from which 
they have been culled ( IK06I) . Apart from that, galaxies in 
our sample are distributed over similar ranges of physical 
properties (e.g., mass, SFR, extinction) compared to other 
samples based on rest-frame optical selections. In contrast, 
BX/BM galaxies tend to be less obscured by dust as ex- 
pected from their selection technique and the fraction of ob- 
jects with strong Ha flux and luminosity (e.g., /A(Ha) > 
10~'^erg s"'cm"^ and L(Ha) > 2 x lO'^^erg s"') appears to be 
smaller than among iBz/Ts. Since dust extinction is not cor- 
rected in Figure [TO] (c) and (d), the difference between the 
rest-frame optically-selected population and rest-frame UV- 
selected one would become even larger once extinction cor- 
rection is applied. 

5.3. Comparison of Star- formation Rates from UV, Ha, and 

Mid-infrared 

In 34.21 we have derived Ha star-formation rates with ex- 
tinction correction following, in which the dust extinction to- 
wards H II regions is larger than that resulting from the SED 
fit to the stellar continuum. To check the validity of this ex- 
tinction correction, SFRs from Ha are compared in Figure [TTI 
with those derived from the MIR and the UV using the same 
procedure as in Daddi et al. (2007a), which is reproduced here 
below. Following the conversion of iKennicu tt ( 1998), the IR 
SFR, or SFR(IR), is derived from the total IR luminosity (Lir) 
as 

SFR(IR) [Moyr-i] = 1.73 x 10-'%r [L©]. (8) 

The total IR luminosity is derived from the luminosity- 
dependent SED Ubrary of lCharv & ElbazI (1200 Ih by using the 
rest-frame 8 /im luminosity (Lg^^m), where the flux density at 
24 fjm from SpitzerlMWSi is used as a proxy of the rest-frame 



8 /im, since 24 /im corresponds to ~ 8 /im at z ~ 2. The 
conversion from Lg^m to Lir is then 

log (g)=1.501og(^) -4.31, (9) 

if log(j/L8^ni) > 9.75, and 

^°Kz^)=^-'^^^°K^)^^'^^' ^^^^ 

if log(i/L8p,n) < 9.75 (iDaddi et al.l l2007ah . The UV SFR, 
or SFR(UV), is derived from rest- frame 1500 A l uminosity 
(Lisoo) by using equation (5) of iDaddi et al.l (l2004al) . i.e., 

SFR(UV) [Mgyr-'l = 1.13 x lO-^^Lisoo [ergs s"!]. (11) 

The observed-frame B-band is used to derive L1500 after ap- 
plying a /T-correction based on the redshift and the UV- 
slope of each object. Finally, the reddening E(B-V) from 
Equation ([T|i is used for the correction of dust extinction 
to derive the extinction-corrected SFR. Figure fTTT a) com- 
pares the total SFRs derived from summing SFR(UV) uncor- 
rected for extinction and SFR(IR), and SFRs derived from 
Ha luminosities with the extinction correction described in 
34.21 These two sets of SFRs agree well with each other, 
with a few exceptions. One sBzK, D3a- 11391, having much 
lower SFR(Ha) = 4OM0yr-' comp ared to SFR(UV-hI R) = 
1 lOOMoyr"' is a MIR-excess object (IDaddi et alJl2007bh . de- 
fined as an object having SFR(UVh-IR) > 3 times higher 
than the SFR(UV) corrected for extinction. The object has 
extinction corrected SFR(UV) of 120M(7)yr~' which i s also 
much smaller than SFR(UVh-IR). IDaddi et alj (l2007bl) sug- 
gested that MIR-excess objects could be hosting an obscured 
AGN, being almost completely opaque in the UV due to dust 
extinction, but emitting strongly in the MIR from hot dust sur- 
rounding the nuclei. 

The outliers showing very large excess in Ha SFR, D3a- 
5814 and D3a-6397, are possibly due to overestimates in 
E{B-V) from the SED fitting. The E{B-V) of these 2 sBzKa 
from BzK colors are smaller than those from SED fitting (Ta- 
ble [T] and |4|i. For large values of E{B-V), even a relatively 
small difference in E(B-V) makes a large difference in the 
resulting Ha SFR. 

In the original Calzetti's recipes for extinction correction 
it is suggested that the emission lines from H II regions suffer 
more extinction than the stellar continuum by a constant factor 
of2.3, or£'(B-y),s,ar = 0.44£'(B-y)ga,s (Calzetti 2001). If we 
simply use £'(Z?-y)stai from SED fitting as £'(Z?-y)gas, i.e., no 
additional obscuration toward H II regions, then the resulting 
Ha SFRs are lower than those derived in 34.21 bv factors of 
- 2 to 5. On the other hand, \fE{B-V\,^r = OAAE{B-V)^^^ 
is used for extinction correction following Calzetti ( 2001), it 
generally overest imates the Ha SFR up to factor of ^ 2. Ha 
SFRs derived by iHayashi elai.. (2009) by using E{B-V)^^^ 
according to the Calzetti's recipes are overestimated by a large 
factor compared to the UV SFRs with extinction correction 
by £■(/?- y)stai. On the other hand, extinction corrected Ha 
SFRs and UV SFRs agree reasonably well for the UV-selected 
z ~ 2 galaxies (Erb et al. 2006c) in which the amount of dust 
extinction is not large on average. This suggests that original 
Calzetti law could produce over-correction at least for heavily 
obscured galaxies like sBzKs. Therefore, the two cases above, 
£(B-V)gas =£(B-y)star and £(B - y)ga,s = £(B " y)star/0.44, 
bracket the Ha SFR derived in 34.21 bv using the recipe of 
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ICid Fernandes et aTl ( 120051) . and additional extinction toward 
H II regions depending on the amount of extinction of stellar 
components could be justified in correcting the Ha fluxes. 

In Figure [TT]UV SFR corrected for dust extinction and Ha 
SFR are plotted. After extinction correction for the Ha lu- 
minosities, both SFRs agree reasonably well except for D3a- 
5814 which is also an outlier in Figure [TTTa). The other out- 
lier in Figure fTTl a). D3a-6397, is not an outlier in Figure 
[TTT c). possibly because of the smaller difference between the 
E{B-V) from SED fitting and that from BzK color, compared 
tothatofD3a-5814. 



5.4. Broad Ha Emission Lines and Super Massive Black 

Holes 

As shown in the above sections, D3a- 11391 is the only 
galaxy in our sample of Ha detected BzK source s that would 
be classified as a MIR-excess galaxy following Daddi etaP 
(!2007b). This galaxy shows a broad Ha component in ad- 
dition to the narrow component, supporting the idea that 
it contains a powerful AGN. The velocity width of D3a- 
11 391, FWHM ~ 2450 km s~', is in between th at of stacked 
SINS galaxies at z ~ 2 ("Shapir o et all l2009h with > 
IO'^Mq (~ 2200 km s"') and "that of stacked SINS AGNs 
(c± 2900 km s ~'). Many of the most massive galaxies in 
Shapiro et al.' ("2009") sample are also MIR excess galaxies. 
Swinbanket al. (2004) also detected a broad-line component 
in the stacked spectrum of submillimeter galaxies (SMGs), 
with FWHM = 1300 km s~' for the stacked spectrum of the 
whole SMG sample, and FWHM = 890 km s"' for fliat of the 
SMGs with no sign of an AGN component. Moreover, some 
individual SMG s show a broad Ha with FWHM c:^ 2000- 
4000 km s-i (Swinbank et"!!] 120041: [Alexander et al.1 [20081) . 
The velocity width of the broad-line component in our stacked 
spectrum, FWHM = 870 km s~'(cf. ^4.5l l is close to the veloc- 
ity width of the SMGs without obvious signs of AGN activity. 

Although strong supernova-driven winds can produce a ve- 
locity width of FWHM > 2000 km s"', it seems difficult to ex- 
plain the MIR-excess as due to supernova remnants because 
supernovae are directly related to the star-formation activity 
which would also produce UV and Ha emissions, consistent 
with the MIR emission. Thus, the object could be hosting 
an AGN at the center and the MIR-excess can be due to the 
surrounding dust torus which absorbs the UV emission from 
the central nucleus. Although typical MIR-excess galaxies 
are expected to contain very obscured, often Compton thick 
AGNs, the broad Ha emission can be explained as a leak 
from the broad-line region (BLR) around the central nucleus, 
which would still be observable in the optical/NIR where the 
optical depths is lower than in the UV. Alternatively, Ha pho- 
tons in the broad component could have been scattered into 
the line-of-sight by either electrons or dust particles. Unfor- 
tunately, there are no X-ray d ata for this object , and t he Ha- 
X-ray diagnostics discussed inl Alexander et al.[ (120081) cannot 
be applied. In the following discussion, we assume that the 
broad-line emission is not dominated by scattered light but is 
seen directly. 

By using the measured Ha luminosity and the width of the 
broad line and assuming that it is coming from an AGN, we 
attempt here a rough estimate of the virial mass of the central 
super massive black hole (SMBH) by adopting the relation 



provided bv [Greene & H0I (l2005h : 



Mbh = {2.0^03) X 10' 



0.55±0.02 



lO^^ergs s- 



/FWHMHa 

VlO^kms-i 



2.()6±0.06 



(12) 



or the relation by iKaspi et al] ( 120051) as converted by 
[Greene & Hoi ( [2005h for the Ha luminosity and width: 



Mbh = (1.3 ±0.3) X 10" 

().57±0.06 



1044ergs s- 



/FWHMna 
VlO^kms-i 



2.06±0.06 



Ma 



(13) 



For the broad component of D3a-11391, Lhq = 7.6 x 10"*^ 
ergs s"' and FWHMna = 2450km s"' lead to Mbh ^ 3 x 
10^ M0 by averaging the values from these two equations. 
If the geometry of BLR is disk-like, the average correction 
factor for the inclination wo uld be ~ 2.7 ([McLure & DunlopI 
[2002; Alex ander et al.l2008h . which gives Mbh ~ 9 x lO^M©. 
The resulting value is ~ 2-3 times smaller than the average 
SMBH mass of broad-line SMGs at z ~ 2 ([Alexander et alj 
[20081) . and a factor > 10 smaller than that of optically se- 
lectedQSOsinSpSSa t z = 1.8-2. 1 ( [McLure & Dunlop 200-1 
[Alexander eTall [2008[) . On the other hand, SINS galax- 
ies with similar stellar mass have ~ 3 times smaller black 
hole mass compared to what estimated here for D3a- 11391 
([Shapiro et alj[2009l) . Dust extinction for the broad-line Ha 
emission is not considered here since H/3 emission line is out 
of the observed wavelength range. If we assume the same 
amount of extinction as that for the host galaxy {Ay = 1), it 
would increase the virial black hole mass by about a factor of 
2. 

Rest-frame 2-10 keV X-ray luminosity can be inferred 
from rest-frame 8 /j,m luminosity assumi ng that the emis- 
sion from AGN dom inates at 8 fj,m ([Lutz et al.[ [lool 
[Alexander et~an I2OO8I) . Since 24 /im corresponds to rest- 
frame 8.7 ^m at z = 1.774, MIPS 24 /im flux can be used 
as a proxy of rest-frame 8/xm luminosity. The flux of D3a- 
1 1391 is ~ 400 /iJy at 24 /im, corresponding to vL^^^ ~ lO'*^ 
ergs s~', which is translated into an absorption corrected 2-10 
keV luminosi ty of L^-micRV — 1 0^"^ er gs s"' by using the same 
method as in [Alexander et al] ([2005i) . Alternatively, the Ha 
luminosity c an be also used to derive 2-10 keV luminosi 
( [Ward et al.[[ 19881) . and from Figure 5 of [Ward et al] ( [198 
we derive L2-iokeV — 10'*'* ergs s"', in agreement with the es- 
timate from the 8 /im luminosity. Given that a mass-accretion 
of 1 Mq yr"' corresponds to Li-iokev — 2.2 x lO"*'* ergs s"' , the 
mass-accretion rate of matter to the central SMBH of D3a- 
11391 can be estimated as ~ O.5M0 yr"'. While this esti- 
mate is obviously affected by large uncertainties, we note that 
this SMBH accretion rate is about 1/2 of th e value estimated 
for broad-line SMGs ([Alexander et alJl2008[) and one order of 
magnitude smaller than the typical value of QSOs a t z= 18- 
2.1 (.McLure & Dunlop[|2004t [Alexander et atll^OS '). 

Then we can proceed further, and crudely estimate the Ed- 
dington fac tor, r?, by comparing Mbh and the accretion rate 
(Figure 2 of [Alexander et al. 2008). For D3a-1 1391 this gives 
rj ~ 0.3-0.9, depending on the assumed geometry, somewhat 
higher than the typical Eddington ratio of broad-line SMGs 
and lower than the average of SINS galaxies. However, un- 
certainties in both black hole mass and mass-accretion rate are 
large. 
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A broad-line component with FWHM = 870 km s"' was 
also detected our stacked spectrum of the remaining BzK 
galaxies (Fig. |7]i, although not with high significance. Its ve- 
locity width is not as large as that of BLRs, but rather similar 
to that of narrow-line Seyfert galaxies. Proceeding in the same 
way as for D3a-1 1391, we then estimate the black hole mass. 
X-ray luminosity, accretion rate and Eddington factor for the 
average galaxy represented by the stacked spectrum. We 
then obtain Mbh = 3 x 10* or 7 x lO^M©, respectively for 
spherical and disk-like geometry, L2_iokeV — 5 x lO^-' ergs s"', 
an accretion rate of Q.2Mq yr"', and ij > I. Here we have 
used the average redshift (z = 1 .96), average Ha luminosity 
(3.4 X lO''^ ergs s"' cm"^) and average 24fim flux (156/iJy) 
of the galaxies in the stack. It is also assumed that the 
whole 24/im flux comes from an obscured AGN, which would 
overestimate the rest-frame 8/im luminosity. By inferring 
the X-ray luminosity from the broad Ha luminosity we get 
Li-iokev — 3 X 10"^^ erg s"', quite consistent with the value 
estimated fr om the 8/im lumino sity, and similar to the value 
estimated by lDaddi et al.l (l20071j) for the MIR-excess galaxies 
on the basis of X-ray stacking. 



5.5. Mass-Metallicity Relation 

The gas-phase oxygen abundances of the galaxies in the 
present sample as derived from Equation (|5]l are plotted as 
a function of stellar mass in Figure [12] Error bars for stel- 
lar mass are set to 0.3 dex, as the typical uncertainty of 
our estimates. AGN candidates, i.e., those showing large 
[N II]/Ha (> 0.7) ratios, are also plotted in Figure [12] as 
open squares and circles. The metallicity derived from the 
stacked spectrum (with and without the broad line compo- 
nent) is also shown along with the average stellar mass of the 
stacked galaxies (star symbols). The solar oxygen abundance 
is indicated with the dashed line, hence almost all galaxies 
in the present sample appear to have super-solar ISM abun- 
dances. However, the absolute values of the oxygen abun- 
dances should be taken with caution having been obtained 
indirectly from the [N II]/Ha ratio. For example, by com- 
paring various indicator s of gas-phase ox ygen abundance, in- 
cluding the T; method, IShi et al.l (120051) using [O III]A4363 
found that in the luminosity-metallicity (L-Z) relation of the 
local blue compact dwarfs, there is a systematic offset be- 
tween the L-Z relation derived from the method and that 
from the 7?23 method, i.e., the method that we have used as 
a reference calibration by scali ng t he ionization parameter in 
the N2 -based calibration (see 34.3l l. The /?23-based L-Z rela- 
tion has a systematically higher zeropoint than the Tg-based 
one, thus the Rxi method may overestimate the true oxygen 
abundance. Therefore, using the same calibrator is crucial to 
compare various observations as shown below. 

,T04. derived the M^-Z relation in the loc al universe (z ^ . 1 ) 
by using 53,000 star-forming galaxies from the SDSS DR2 
release, and their M*-Z relation can be used as the refer- 
ence relation for z — Q- However, their metallicity measure- 
ments are based on model fitting to multiple emission lines 
(jCharlot & Longhetti 2001), and therefore are not directly 
comparable to our results derived from Equation ([5]). To cope 
wit h thi s problem, we have used SDSS DR4 data as described 
in 34.31 and multiplied their stellar masses by a factor 1.5 (to 
convert from Kroupa IMF to Salpeter IMF). Then the result- 
ing M^-Z relation for SDSS galaxies with Rn 10^,2 calibrated 



abundances is fitted with the 2nd order polynomial, 

12 + log(0/H)= 1.0512 

+ 1.38361og(M,/M0) (14) 

-0.0602 [log(M,/M0)]', 

which is finally used as the reference relation at z ^ in Fig- 
ure[l3] 

The My,-Z relation for 57 star-forming galaxies from the 
GDDS/CFRS surveys (Savaglio et al. 2005) is then used 
as a reference at intermediate redshift, z — 0.7. Since 
they used Rn as in iKK04. to calibrate abundances and the 
iBaldrv & Glazebrookl ( l2003h IMF, we convert only their stel- 
lar masses to those corresponding to the Salpeter IMF, by mul- 
tiplying them by a factor 1.8. 

In addition to the sBzKs in this study, we also consider 
masses and abundances of the following z ^ 2 objects: 6 dis- 
tant r ed galaxies (DRGs) at 2.4 < z < 3.2 (van Dokkum et aj] 
|2004|) where only one out of 6 presented in lvan Dokkum et al] 
(2004) is used here because 4 DRGs show obvious AGN 
features (broad-line and high [N II]/Ha ratio), and for an- 
other object neither N2 nor R23 indicators could be used 
for a metallicity measurement; 87 UV-selected (B X/BM) 
star-forming galaxies at (z) = 2.26 dErb et alj|2006al) which 
are binned into 6 mass bins containing ~ 15 galaxies per 
bin; the Lyman-break galaxy MS1512-cB58 at z = 2.73 
(Teplitz et al. 2000; Baker etal. 2004); and final ly die aver- 
age of 7 SMGs with [N II]/Ha < 0.7 at (z) = 2.4 ( Small et alj 
i2004j; iSwinbank et al, 20 04). The stellar masses are all cor- 
rected to the Salpeter IMF, taking gravitational lensing into 
account in the case of cB58. All gas-phase oxygen abun- 
dances are derived with the KK04 calibration by using the 
[N II]/Ha emission line ratios as listed in the quoted refer- 
ences so to ensure full homogeneity with the sBzK abun- 
dances derived in this paper. 

Figure [T3] compares masses and abundances of sBzKs at 
z ~ 2 to those of the objects mentioned above. Most sBzKs 
at z — 2 are already enriched to abundances comparable to 
those of SDSS and GDDS/CFRS galaxies in the same stel- 
lar mass range, i.e., logM^, > 10.5. The ISM abundances of 
sBzKs appears to be consistent with those of z — 2 galaxies 
of comparable mass, although selected according to differ- 
ent criteria (i.e., DRGs and S MGs), while the M^-Z relation 
of UV-selected galaxies from lErb et alj (l2006al) tend to show 
slightly lower (~ 0.15-0.2 dex) abundances at a given stel- 
lar mass. This offset from the M*-Z relation of UV-selected 
galaxies which is generated by the spectral stacking can be 
seen also when compared with the metallicities of the com- 
posite spectra of sBzKs at similar stellar mass. Since the 
metallicity uncertainties are rather large, this offset may not 
be very significant. Also Havashi et al. (2009) have reported 
a ~ 0.2 dex higher average metallicities of their sBzKs com- 
pared to U V-selected galaxies. Our metallicities are closer 
to flios e of lHavashi et alj (120091) than to diose of lErbetalJ 
(l2006al) . though uncertainties are large. iHavashi et all (l2009t) 
also mentioned a possible bias towards higher metallicity ob- 
jects, which would enable us to detect [N II] emissions, hence 
to measure the metallicity. Indeed, their stacked spectrum of 
[N II] undetected objects yields roughly the same metallicity 
a s that of UV-selected g alaxies. 

ISavaglio et alj (120051) noted the differential redshift evolu- 
tion of the M+-Z relation from z^0.7toz = 0.1, with the most 
massive galaxies at z ~ 0.7 (logM^ > 10.3) being already en- 
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riched to the ISM abundances of the local galaxies, while less 
massive galaxies start leaving the local M*-Z relation already 
at redshift as low as z ~ 0.1. Our result extends this trend to 
higher redshifts, finding that at the high mass end the chemi- 
cal enrichment was virtually complete by z '--^ 2. This indicates 
that massive galaxies evolve faster than less massive galaxies, 
where star formation and ensuing chemical enrichment are 
slow and last longer, i.e., yet another manifestation of down- 
sizing. 

5.6. Relation between Stellar Masses and Specific 
Star-formation Rates 

The correlation between SFR and stellar mass has been 
extensively in vestigated from the nearby universe to high 



redshift (e.g.. iBrinchmann et al 
Elbaz et al.l 120071; iNoeske et al 



2004 



Pannella et al.l 120091; iSantini et al. 



20071 
120091) 



IDaddi et aLipOOTa : 
iDunne etal.1 12009 : 
In particular, the 
specific star-formation rate (SSFR), defined as the SFR per 
unit stellar mass, is widely used to quantify the contribution 
of current star-formation activity to the growth of the total 
stellar mass, i.e., of how efficiently stars are formed. 

In Figure [14] we plot the SSFR as a function of stel- 
lar mass for sBzK galaxies with Ha detection. Though the 
scatter and uncertainty in stellar mass are large, the average 
SSFR (in yr"') is around log SSFR :i; -9, and the SSFR-M,^ 
does not show a detectable slope. A comparison of M*- 
SSFR relations from different sources and at different red- 
shifts is also shown in Figure [15] Data from literature are 
taken from Elb az et al.l (12007 ) for star-forming galaxies at 
z ~ in SDSS and z ~ 1 in GOOD S, as well as z ~ 2 
objects including sBzKs (Da ddi et al.l (2007a; Pannella et id] 
[20091) . UV-selected star-forming galaxies (.Erb et al...2006dK 
DRGs dya n Dokkum et al. 2004), and average of SMGs 
dSmail etal. 2004; Swinb ank et al.li2004 . lElbaz et all dlOOTI) 
and Daddi et al. (2007a) have estimated SFRs by adding ex- 
tinction uncorrected rest-frame UV and IR luminosities to 
trace both absorbed and unabsorbed star fo rmation. SFRs for 
sBzKs in the study of lPannella et al.l (120091) have been derived 
from 1.4 GHz radio data (hence independent of extinction) 
finding excellent agreement with Daddi et al. (2007a). SFRs 
of DRGs and UV-selected star-forming galaxies at z — 2 are 
derived fr om Ha luminosities with e xtinc tion corre c tion de - 
scribed in Ivan Dokkum et all (120041) and lErb et al.1 (l2006cl) . 
respectively. The SFR of SMGs is derived from sub-mm 
emission. Since these SFRs are taken from the original pa- 
pers, which besides different selection criteria have also used 
a variety of different observables (such as UV-luminosity, Ha 
luminosity, FIR luminosity and radio luminosity), it is not a 
surprise to find even large systematic differences. 

From low to high redshift, SSFRs at a given stell a r mas s 
increase systematically as shown by e.g.. El baz et al.l (120071) . 
IDaddi e t al. (2007a) and Pannella et al. (200l). The sBzKs 
presented in this study distribute around the average SSFR- 
relation for sBzKs in these two latter studies. One sBzK 
galaxy shows a SSFR as high as that typical of SMGs, and can 
be regarded as a real outlier, though it is also possible that its 
Ha SFR may have been overestimated since the object shows 
an Ha excess but appears normal in the SFR(UV)-SFR(IR) 
comparison (Figure [TTTi. Note that the M*-SSFR relations for 
SDSS star-forming galaxies at z 0, z — 1 and for sBzKs at 
z ~ 2 have almost flat slopes, ranging from ~0 to -0.23, while 
the distribution of UV-selected galaxies in lErb et al.] (I2006ah 
shows a much steeper slope. A steeper slope, but systemat- 
ically higher SSFRs for sBzKs in the M^-SSFR relation has 



been also reported by iHavashi et al.l ( 1200 91). having used Ha 
SFRs. However, their Ha SFR are systematically higher than 
SFRs from the UV, while no co mpar ison with the total (e.g., 
UVh-IR) SFR is available (see and ^53] l. On average^^ 
sBzKs in the present stud y, as wel l as in those of Dad di et al] 
(2007a), Pannell a et al.. (.20091) and lHavashi e t al. (2 009) , have 
higher SSFRs at the massive end compared to those for UV- 
selected galaxies. 

An almost flat M*-SSFR relation for z ~ 2 Bz^T-selected 
galaxies which is almost indistinguishable from that of 
Daddi et al. (2007a) shown in Fij;ure^[T5] is indeed found by 
Dunne et al.1 (120091) and iPannella et alT t2009). The discrep- 



ancy with respect to other , steep M^-S SFR relations is as- 
cribed to selection effects (iDunne et al.l 2009). e.g., the UV 
selection being biased against massive, highly obscured and 
intensively star-forming galaxies, and in addition to a system- 
atic und erestimate of the dust extinction with increasing mass 
(IPann ella et a l. 2009). Regarding the discrepancy of slopes 
between Ha spectroscopic surveys, an uncertain extinction 
correction for Ha from SED fitting might be one of the ma- 
jor causes. Larger spectroscopic samples, exploring a wider 
range of stellar masses and SFRs, and detecting of H/3 in addi- 
tion to the Ha (so to estimate accurate dust extinctions from 
the Balmer decrement) would be crucial for an independent 
evaluation of the slope in M*-SSFR relation at z — 2. Our re- 
sults, th ough relative to a v ery s mall sample, are in li ne with 
those of IDunne et all (120091) and IPannella et all ( I2009I) . 

5.7. An Interpretation of the Mass-metallicity and 
Mass-SSFR Relations with Simple Evolutionary 
Population Synthesis Models 

In this section we try to interpret the M*-Z and M*-SSFR 
relations with evolutionary population synthesis models. We 
use the PEGASE.2 models that are also used for SED fitting 
to derive stellar masses and dust extinction parameters. These 
are one-zone models which neglect interactions and merging 
of galaxies, that are naturally expected in hierarchical struc- 
ture formation scenario. They also neglect galactic winds 
driven by supernovae and/or AGN feedback, that are known 
to exist and are responsible for the metal enrichment of the 
intergalactic and intracluster media. Although a galaxy rarely 
evolves as is modeled by PEGASE.2, its simple description 
may give us a starting point towards understanding the M*- 
Z and M,t-SSFR relations, and indicate in which direction we 
should try to find more satisfactory solutions. A detailed com- 
parison with the full numerical simulation is not the focus of 
this paper 

Two scenarios, the simple closed-box model and the in- 
fall model for chemical enrichment, are considered here with 
a Salpeter IMF from O.IMr to IIQM.,, and the B-series of 
chemical yields from Wooslev & Weaver ( 1995). These mod- 
els assume that a galaxy is initially a purely gas cloud with 
zero metallicity, Z = 0, and the SFR is assumed to be propor- 
tional to the gas mass (i.e., a Schmidt law with n = 1): 



SFR(f) = — Mg,,(f), 



(15) 



where r^f is the star-formation timescale. Closed-box models 
assume no gas infall and all gas resides in the system at the 
beginning. This is the classical monolithic collapse scenario. 
In the infall models, initially all gas is assumed to be in an 
outer reservoir and falls into the inner star-forming region at 
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a rate of 

Cinfaii (0 = exp ( — ) , (16) 

Tinfall V 'Hnfall/ 

where Tinfaii is the infall time scale. We use models with 
Tsf = 0.1 and 5 Gyr and we assume Tst = rinfaii for infall models. 
Short time scale models with Tsf = 0.1 Gyr can successfully 
reproduce the color-magnitude (C-M) r elations of elliptical 
galax ies in nearby clusters of galaxies jKodama & Arimotol 
119971) . and long time scale models with Tsf = 5 Gyr can explain 
the photometric properties of nearby late type, Sb-Sc, galax- 
ies (Arimoto et al. 1992: iFioc & Rocca-Volmerangd 119971) . 
Following KK04, the ISM oxygen abundance is derived from 
the gas-phase metal mass fraction using the relation: 

12 + log(0/H) = 12 + log(^^^, (17) 

assuming the so lar abundance pattern reported by 

I Anders & G reyess^ (|198 9J and the solar oxygen abun- 
dance from Allende Priet o et al.l (l2001h . 

Figures [16] and [IT] show that both the closed-box and the 
infall models with Tsf = 0. 1 Gyr can roughly reproduce the lo- 
cations of sBzKs and other populations of z ~ 2 galaxies at the 
massive end of the M^,-Z relation and of the Mj,-SSFR relation, 
with 100 < f < 500 Myr ages and 10" < M,otai/M0 < 10'^ 
of gas mass. This range of ages is also roughly consistent 
with t he onset epoch of galactic winds needed by the mod- 
els of IKodama & Arimotol (Il997h to reproduce the C-M re- 
lations of elliptical galaxies in the Virgo and Coma clusters. 
On the other hand, models with longer time scale, e.g., Tsf = 5 
Gyr, cannot fit both relations simultaneously. For example, 
Tsf = 5 Gyr infall models take > 5 Gyr to reach the observed 
metallicities of z ~ 2 galaxies, while the age of the universe 
at z ~ 2 is only 3.4 Gyr for the adopted cosmology. Only 
simple models can reproduce the distribution of sBzKs and 
z ~ 2 galaxies in the M^-Z relation with such long Tsf. How- 
ever, in the Mj,-SSFR relation, it turns out that simple mod- 
els with longer time scale require huge amount of initial gas, 
Motai > 1O'^M0, which is larger than typical stellar mass of 
the most massive early-type galaxies found in the local uni- 
verse. Moreover, the ages and stellar masses mentioned above 
imply an average SFR > lO^Moyr"', much in excess of the 
observed SFRs of galaxies in such mass range. 

On the other hand. Equations ( fTSI ) and ( 1161 ) imply a sec- 
ularly declining SFR, such that all galaxies are assumed to 
have started with their maximum SFR and to be caught at 
their minimum SFR. This assumption is especially doubtful at 
z ~ 2, where the SFR of many star-forming ga laxies may ac- 
tually increase with t ime, rather than decrease (lRenzinill2009l : 
iMaraston et"ani2010h . 

The high SFRs in the BzK galaxies presented here require 
galaxies to host a large amount of molecular gas to be used for 
star formation. Recent millimeter and radio observations have 
been revealing that star-forming BzK galaxies indeed harbor 
larg e amount of CO molecule s, hence also molecular hydro- 
gen (iDaddi et al . 2008, 2010; D annerbauer et alJ2009l) . These 
CO observ ations as well as multi-wavelength SED analysis 
jPaddi et a l. 2007a b) suggest that the gas consumption or 
star-formation timescale could be up to 1 Gyr rather than 
100 Myr inferred in this study. In the framework of PE- 
GASE.2 used here, t < 1 Gyr is in fact an upper limit in or- 
der to reproduce the observed M*-Z and Mj,-SSFR relations at 
the same time, still being consistent with .Daddi et al. (.2007al 



l2008h . 

Again, note that closed-box or infall models are certainly 
oversimplifications over the real star-formation histories of 
galaxies. Recent numerical simulations favor continuous (al- 
beit fluctuating) cold-stream accretion for the main driver 
of the galaxy gro wth, rather than short-lived starbursts (e.g., 
iDekel etaTI 120091) . Thus, infall models may capture this as- 
pect, but the results must critically depend on the assumed 
evolution with time of the infall rate. 

SMGs, which are outlier in the M*-SSFR relation, require 
extremely young ages, a few 10 Myr, while in the M,t-Z rela- 
tion they occupy nearly the same position as other z — 2 star- 
forming galaxies. This is possibly due to vigorous, merger- 
driven starbursts rather than due to internal star-formation 
mode assumed in the one-zone PEGASE.2 models. These 
violent star- formation process es could push a galaxy almost 
vertically (Feulner et al. 2005) in the M*-SSFR diagram, with 
little change of its position in the M*-Z diagram. 

6. SUMMARY AND CONCLUSIONS 

We have conducted NIR spectroscopic observations of 28 
sBzK galaxies at ~ 2 and detected Ha emissions from 14 of 
them. By using the Ha and [N II] lines, we have derived Ha 
SFRs and gas phase oxygen abundances. Stellar masses and 
reddening parameters have been also derived by SED fitting to 
the multi-wavelength photometric data. Our results are sum- 
marized as follows: 

• Stellar masses and reddening parameters derived from 
SED fitting agree we ll with those derived by Bz^^-based 
recipe introduced bv lDaddi et alj (l2004ai) . 

• A comparison of SFRs from different indicators (Ha, 
extinction corrected UV, and UVh-MIR) indicates that 
additional extinction towards H II regions over that de- 
rived from the SED fitting of the stellar continuum 
is required for the Ha SFR to recover the total SFR 
inferred from UVh-MIR. The required additional ex- 
tinction is in agreement with the recipe proposed by 
ICid Fernandes et alJ (120051) . 

• One object, D3a-1 1391, shows MIR-excess and a broad 
Ha component with FWHM ~ 2500km s"' , as well as a 
narrow line component, suggesting that it may host an 
AGN, which would be highly obscured at UV and X- 
ray wavelengths. Although with large uncertainty, we 
have estimated the mass of the SMB H and its intrin - 
sic X-ray luminosity following Alexander et al.l (120081) . 
obtaining Mbh — (3-9) x lO^M©, depending on the as- 
sumed geometry of the accretion disk, and L2-iokeV — 
lO**"* ergs s"'. From these quantities, the mass accre- 
tion rate onto the central SMBH and Eddington fac- 
tor are calculated as Mbh — O.5M0yr"' and 77 ~ 0.3- 
0.9, again depending on the geometry. These proper- 
ties are in between those of normal star-forming galax- 
ies and those of broad-line SMGs, indicating a possible 
evolutionary connection between normal UV/optical- 
selected galaxies, broad-line/MlR-excess sBzKs, and 
star-forming broad-line SMGs that may result from an 
accelerated growth of SMBH. 

• Most of sBzKs presented here have already reached a 
metallicity similar to that of local star-forming galax- 
ies of the same mass. They tend to have higher metal- 
licity compared to UV-selected z ~ 2 galaxies, even at 
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the top mass end. This indicates that sBzKs axe on av- 
erage a more evolved population compared to that of 
UV-selected galaxies at the same redshift. 

• Specific SFRs of most sBzKs are consistent with 
a tight M ^-SSFR coiTelation (Daddietal. 2007a; 
iPannellae t al. 2009), with a couple of outliers reaching 
very high SSFRs, similar to those of SMGs 

• Within the framework of PEGASE.2 closed-box or in- 
fall models, the M^-Z and M*-SSFR relations of z ~ 
2 sBzKs can be reproduced simultaneously assuming 
a very short star-formation or infall timescale, r 2± 
100 Myr, and large gas mass at the beginning. However, 
the implied SFRs averaged over the life of the galaxies 
appears to be excessively large, suggesting that secu- 
larly declining SFRs may not be appropriate for star- 
forming galaxies at z ~ 2. 



We thank the anonymous referee for providing the use- 
ful and constructive report Part of the materials presented in 
this paper is from the Ph.D. thesis of MO at the Astronomy 
Department of the University of Tokyo. We thank Kentaro 
Motohara for providing the reduction pipeline used for the 
OHS/CISCO spectra and for giving advice about data reduc- 
tions. Marcella Brusa is thanked for providing us upper limits 
for the XMM-Newton/Chandra fluxes and Ranga-Ram Chary 
is thanked for help with the MIPS 24 ^m data. We thank the 
staffs of the Subaru telescope and of the VLT for support- 
ing the observations. This work was supported in part by a 
Grant-in-Aid for Science Research (No. 19540245) by the 
Japanese Ministry of Education, Culture, Sports, Science and 
Technology. DMA thanks the Royal Society and Leverhulme 
Trust for financial support. AR express gratitude to the INAF 
- Osservatorio Astronomico di Bologna for its hospitality and 
support. 



REFERENCES 



Alexander, D. M., Smail, I., Bauer, F. E., Chapman, S. C, Blain, A. W., 

Brandt, W. N., & Ivison, R. J. 2005, Nature, 434, 738 
Alexander, D. M., et al. 2008, AJ, 135, 1968 

Allende Prieto, C, Lambert, D. L., & Asplund, M. 2001, ApJ, 556, L63 
Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Arimoto, N., Yoshii, Y, & Takahara, E 1992, A&A, 253, 21 
Baker, A. J., Tacconi, L. J., Genzel, R., Lehnert, M. D., & Lutz, D. 2004, 

ApJ, 604, 125 
Baldry, I. K., & Glazebrook, K. 2003, ApJ, 593, 258 
Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 
Binney, J., & Tremaine, S. 1987, Galactic Dynamics (Princeton Univ. Press) 
Bonnet, H., et al. 2004, ESQ Messenger, 117, 17 

Brinchmann, J., Chariot, S., White, S. D. M., Tremonti, C, Kauffmann, G., 

Heckman. T., & Brinkmann, J. 2004, MNRAS, 351, 1151 
Buat, v., BoselU, A., Gavazzi, G., & Bonfanti, C. 2002, A&A, 383, 801 
Calzetti, D. 2001, PASP, 113, 1449 
Chabrier. G. 2003, PASP, 115, 763 
Chariot, S., & Longhetti, M. 2001. MNRAS, 323, 887 
Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 

Cid Fernandes, R., Mateus, A., Sodre, L., Stasihska, G., & Gomes, J. M. 

2005, MNRAS, 358, 363 
Cimatti, A., et al. 2008, A&A, 482, 21 

Daddi, E., Cimatti, A., Renzini, A., Fontana, A., Mignoli, M., Pozzetti, L., 

Tozzi, P, & Zamorani, G. 2004a, ApJ, 617, 746 
Daddi, E., et al. 2004b, ApJ, 600, LI 27 
Daddi, E., et al. 2005, ApJ, 631, L13 
Daddi, E., et al. 2007a, ApJ, 670, 156 
Daddi, E., et al. 2007b, ApJ, 670, 173 

Daddi, E., Dannerbauer, H., Elbaz, D., Dickinson, M., Monison, G., Stern, 

D., & Ravindranath, S. 2008, ApJ, 673, L21 
Daddi, E., et al. 2010, ApJ, 713, 686 
Dannerbauer, H., et al. 2006, ApJ, 637, L5 



Dannerbauer, H., Daddi, E., Riechers, D. A., Walter, F, Carilli, C. L., 
Dickinson, M., Elbaz, D., & Morrison, G. E. 2009, ApJ, 698, L178 
Dekel, A., et al. 2009, Nature, 457, 451 
Drory, N., Bender, R., & Hopp, U. 2004, ApJ, 616, L103 
Dunne, L., et al. 2009, MNRAS, 394, 3 
Eisenhauer, F, et al. 2003, Proc. SPIE, 4841, 1548 
Elbaz, D., et al. 2007, A&A, 468, 33 

Ellison, S. 2006, The Fabulous Destiny of Galaxies: Bridging Past and 
Present, 53 

Erb, D. K., Shapley, A. E., Pettini, M., Steidel, C. C, Reddy, N. A., & 

Adelberger, K. L. 2006a, ApJ, 644, 813 
Erb, D. K., Steidel, C. C, Shapley, A. E., Pettini, M., Reddy, N. A., & 

Adelberger, K. L. 2006b, ApJ, 646, 107 
Erb, D. K., Steidel, C. C, Shapley, A. E., Pettini, M., Reddy, N. A., & 

Adelberger, K. L. 2006c, ApJ, 647, 128 
Feulner, G., Gabasch, A., Salvato, M., Drory, N., Hopp, U., & Bender, R. 

2005, ApJ, 633, L9 
Fioc, M., & Rocca-Volmerange, B. 1997, A&A, 326, 950 
Fioc, M., & Rocca-Volmerange, B. 1999, astro-ph/99 12179 
Forster Schreiber, N. M., et al. 2009, ApJ, /Ub, 1 Jb4 
Genzel, R., et al. 2006, Nature, 442, 786 
Greene, J. E., & Ho, L. C. 2005, ApJ, 630, 122 
Hayashi, M., et al. 2009, ApJ, 691, 140 

Iwamuro, F, Motohara, K., Maihara, T., Hata. R., & Harashima, T. 2001, 

PASJ 53 355 
Kajisaw'a, M., et al. 2009, ApJ, 702, 1393 

Kaspi, S., Maoz, D., Netzer, H., Peterson, B. M., Vestergaard, M.. & 

Jannuzi, B. T. 2005, ApJ, 629, 61 
Kauffmann, G., et al. 2003, MNRAS, 346, 1055 
Kennicutt, R. C. 1998, ARA&A, 36, 189 
Kewley, L. J., & Dopita. M. A. 2002, ApJS, 142, 35 
Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183 
Kobulnicky, H. A., & Kewley, L. J. 2004, ApJ, 617, 240 (KK04) 
Kodama, T., & Aiimoto, N. 1997, A&A, 320, 41 
Kong, X., et al. 2006, ApJ, 638,, 72 (K06) 
Kriek, M., et al. 2006a, ApJ, 645, 44 
Kriek, M., et al. 2006b, ApJ, 649, L71 
Kriek, M., et al. 2007, ApJ, 669, 776 

Kriek, M., van der Wei, A., van Dokkum, P. G., Franx, M., & lUingworth, 

G. D. 2008a, ApJ, 682, 896 
Kriek, M., et al. 2008b, ApJ, 677, 219 
Kroupa, P 2002, Science, 295, 82 
Lanzoni, B., & Ciotti, L. 2003, A&A, 404, 819 
Le Fevre, O., et al. 2005, A&A, 439, 845 
Lilly, S. J., et al. 2007, ApJS, 172, 70 

Lutz, D., Maiolino, R., Spoon, H. W. W.. & Moorwood, A. F. M. 2004, 

A&A, 418, 465 
Madau, R 1995, ApJ, 441, 18 

Madau, P., Ferguson, H. C, Dickinson, M. E., Giavalisco, M., Steidel, C. C, 

& Fruchter, A. 1996, MNRAS, 283, 1388 
Magdis, G. E., Rigopoulou, D., Huang, J.-S., & Fazio, G. G. 2009, MNRAS, 

1768 

Maraston, C. 1998, MNRAS, 300, 872 

Maraston, C. 2005, MNRAS, 362, 799 

Maraston, C. et al. 2010, ApJ, submitted 

McCracken, H. J., et al. 2010, ApJ, 708, 202 

McLure, R. J., & Dunlop, J. S. 2002, MNRAS, 331, 795 

McLure, R. J., & Dunlop, J. S. 2004, MNRAS, 352, 1390 

Mignoh, M., et al. 2005, A&A, 437. 883 

Motohara, K., et al. 2001, PASJ, 53, 459 

Motohara, K., et al. 2002, PASJ, 54, 315 

Noeske, K. G., et al. 2007, ApJ, 660, L43 

Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 

Osterbrock, D. E., & Ferland, G. J. 2005, Astrophysics of Gaseous Nebulae 
and Active Galactic Nuclei (2nd ed.; Univ. Science Books) 

Pagel, B. E. J., Edmunds, M. G., Blackwell, D. E., Chun, M. S., & Smith, G. 
1979, MNRAS, 189, 95 

Pannella, M., et al. 2009, ApJ, 698, LI 16 

Pettini, M., Shapley, A. E., Steidel, C. C, Cuby, J.-G., Dickinson, M., 
Moorwood, A. F. M., Adelberger, K. L., & Giavalisco, M. 2001, ApJ, 
554, 981 

Popesso, R, et al. 2009, A&A, 494, 443 

Reddy, N. A., Erb, D. K., Steidel, C. C, Shapley, A. E., Adelberger, K. L., & 

Pettini, M. 2005, ApJ, 633, 748 
Renzini, A. 2005, The Initial Mass Function 50 Years Later, 327, 221 
Renzini, A. 2009, MNRAS, 398, L58 
Renzini, A., & Daddi, E. 2009, The Messenger, 137, 41 
Salpeter, E. E. 1955, ApJ, 121, 161 
Santini, P, et al. 2009, A&A, 504, 751 
Savaglio, S., et al. 2005, ApJ, 635, 260 
Shapiro, K. L., et al. 2009, ApJ, 701, 955 
Shi, F, Kong, X., Li, C, & Cheng, F Z. 2005. A&A, 437, 849 
Smail, L, Chapman, S. C, Blain, A. W., & Ivison, R. J. 2004, ApJ, 616, 71 
Steidel, C. C, Shapley, A. E., Pettini, M., Adelberger, K. L., Erb, D. K., 

Reddy, N. A., & Hunt, M. P 2004, ApJ, 604, 534 
Storchi-Bergmaim, T., Calzetti, D., & Kinney, A. L. 1994, ApJ, 429, 572 



NIR Spectroscopy of sBzKs at z ~ 2 



13 



Swinbank, A. M., Smail, I., Chapman, S. C, Blain, A. W., Ivison, R. J., & 

Keel, W. C. 2004, ApJ, 617, 64 
Teplitz, H. I., et al. 2000, ApJ, 533, L65 
Tremonti, C. A., et al. 2004, ApJ, 613, 898 
van Dokkum, P. G., et al. 2004, ApJ, 611, 703 
Vanzella, E., et al. 2005, A&A, 434, 53 
Veilleux, S., & Osterbrock, D. E. 1987, ApJS, 63, 295 



Ward, M. J., Done, C, Fabian, A. C, Tennant, A. P., & Shafer, R. A. 1988, 

ApJ, 324, 767 
Woosley, S. E., & Weaver, T. A. 1995, ApJS, 101, 181 
Yuan, T.-T., & Kewley, L. J. 2009, ApJ, 699, L161 



14 Onodera et al. 

D3a3287 D3a4626 D3a4654 D3a4751 D3a5814 D3a6397 D3a7429 



i 

# 

■ « 




\ 




■ 


* 


■ 




ililtf^ii •^rriT 




1 ■'V'i''i'ritfi ra'i^i 




■ '"t ■ ■■ 

■ -m 







D3a8608 D3a11391 D3a12556 



Dad1901 Dad2426 Dad3551 



ff 


■ 


p 






IBS' ■ ■* 




Figure 1. Cutout images (~ 15" X 15") of sBzKs with Ho detection in our sample. North is up and east is to the left. Each pair consists of composite color 
image (top) and AT-band image (bottom). Composite color images are created with Subaru/Suprime-Cam BRz' and Biz' for Deep3a and Daddi fields, respectively. 
Objects are located at the center except for Dad 2426-b which can be seen at 3 arcsec west from Dad 2426. 
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Figure 2. One-dimensional spectra around the Ha emission line for the spectra taken with OHS/CISCO. Black, orange, and green lines correspond to raw, 40 
A boxcar smoothed, and Icr sky fluctuation spectra, respectively. Gaussian functions fitted to the smoothed spectrum ai'e also shown with dashed lines which 
respectively correspond to Ha (light green), [N II]A6548 (magenta) and [N II]A6583 (light blue). The sum of these three Gaussian is shown with blue solid line. 
Positions of [N II]AA6548,6583 and Ha are indicated by gray dashed fines with circles and squares, respectively. 
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Figure 4. The SEDs of Ha detected sBzKs: filled squares siiow the observed flux from broad-band photometry, while solid lines show the best-fit SED from 
PEGASE.2 models. Open circles show model magnitudes derived by convolving best-fit SED with the filter response curves. The agreement between the 
observed and fitted photometry is typically very good and so often the open circles are hidden by the filled circles. Note that Dad 2426-b is excluded here because 
it is a serendipitous object not in our iT-selected catalog. 
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Figure 6. The correlation between the metalhcity derived from the R23 / O32 parameter (Equation|6) and from the N2 indices (Equation|5j for 75,561 star-forming 
galaxies in the SDSS DR4 release. The various panels refer to different assumed values for the ionization parameter q in the N2 calibration, as indicated. It 
appears that ~ 3 X lO' is appropriate for metal-rich galaxies if a constant value is assumed. 
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Figure 7. The stacked spectram of the six SINFONI Hct-detected objects (soUd black Unes). The positions of Hq and the two [N II] hnes are marked at the 
top of the panel with dashed lines. Best fit multi-Gaussian profiles are shown with solid blue lines for the total spectrum, and with solid magenta lines for its 
individual components. In the top panel, a broad Ha component is included while only the narrow line components are used for the fit in the bottom panel. 
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Figure 8. The {B-z) vs. (z-A") colors of the objects in the Deep3a field and the Daddi field. Black dots show objects in the A'-selected catalog, while open circles 
(red) and open triangles (blue) with sizes proportional to the X^-band luminosity show sBzKs in the present study, with and without Ha detection, respectively. 
Solid, dotted, and dashed lines separate sBzKs, passive BzK galaxies (pBzKs), and stars (below the dashed line) from other populations (see iDaddi et al.H2004bl : 
IKong et aU200g) . 
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Figure 9. Top: the B^A'-based E(B — V) as a function of BzK-based stellar mass. Symbols are as in Figure[8] Bottom: the same as in the top panel, but for the 
extinction con'ected SFR derived from the rest-frame 1500A luminosity, with extinction correction based on the E(B-V) values plotted on the top panel. 
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Figure 10. Stellar mass, reddening and Ha properties as a function of the X'-band magnitude are plotted for the galaxies in the present study and those from 
other Ha surveys. The various panels give: (a) stellar mass, (b) reddening, (c) Ha flux, and (d) Ha luminosity. Our sBzK sample is shown with filled indigo 
squares, rest-frame UV-selected BX/BM galaxies from .Erb et al j Ii2006b c) and rest-frame optically selected sISzKs from Hayashi et al. 1 2 00^) are shown with 
pin k triangles and gray upside-down triangles, respectively. Circles represent z ~ 2 star-forming galaxies observed by SINS survey I Fo rster Schreiber et al.l 
I200f^) . For SINS galaxies. BX/BM galaxies and rest-frame optically selected galaxies are shown with pink and gray symbols, respectively. 
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Figure 11. Comparisons of SFRs from different indicators: (a) the SFR from Hq vs. the SFR from extinction uncorrected rest-frame UV luminosity and 
far-infrared luminosity deiived from the 24/xm flux: (b) the SFR from rest-frame UV luminosity vs. SFR (UV-l-IR); (c) the SFR(Ha) vs. extinction-corrected 
SFR(UV). In all panels, open and filled symbols refer to SFRs before and after extinction correction, respectively, connected by vertical dotted lines. An object 
with MIR-excess is shown with square symbol and objects with SFR(Ha) much higher than other estimators are shown as triangles. 
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Figure 12. The M^-Z relation for sBzKs at z ~ 2. Filled and open circles and squares represent objects with [N Ilj/Ho < 0.7 and [N II]/Ha > 0.7, respectively. 
Circles and squares represent metalHcities from OHS/CISCO and SINFONI data, respectively. A typical eiTor bar for the stellar mass of 0.3 dex is indicated for 
each galaxy. Stars represent the values derived from the stacked spectra by means of profile fitting with a nd without a broad-line Ha component (open and filled 
star, respectively). The horizontal dashed line shows the solar oxygen abundance, 124-log(0/H)Q = 8.69 4AIlende Prieto et aU200H) . 



24 



Onodera et al. 




8 9 10 11 12 

log(M^/M®) 

Figure 13. The M*-Z relations for galaxies at z = 0.1, ~ 0.7, and ~ 2 from various samples and sources. Filled indigo squares with border show the non- 
AGN sBzKs in our sample observed by OHS/CISCO and SINFONI, and filled star and the open stars refer to the stacked spectrum, as in Figure [T2I Solid and 
dotted orange lines show the local M»-Z relation derived from SDSS data (see text), along with its ilcr range. Green filled triangles represent z — 0.7 galaxies in 
GDDS/CFRS sample I Savaglio et al. 2005). A DRG ( van Dokkum et al. 2004) is plotted as a magenta upside-down triangle. The mean Mt-Z relation for BX/BM 
galaxies is shown by the pink dot-dashed line. Cyan circle and red square represent MS 1512-cB58 (Teplitz et al. 2000; Baker et al. 2004) and the median value 
for the SMGs cSmail et alj,2004. : .Swinbank et al.ii2004) , respectively. In this plot eiTor bars are omitted and stellar mass and metallicity are converted to the 
Salpeter IMF and to the lKK04l rnetallicitv calibration, respectively. The horizontal dashed line indicates the solar oxygen abundance. 
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Figure 14. The SSFRs of sBzKs in this study are plotted as a function of galaxy stellar mass. Symbols are the same as in Figure [T2I A typical mass error of a 
factor of 2 (0.3 dex) is adopted for all galaxies. The objects possibly dominated by an AGN are not plotted in this figure. 
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Figure 15. The SSFRs for different galaxy samples are plotted as a function of galaxy stellar m ass. Filled i ndigo squares with border show non-AGN sBzKs 
in our sample. Magenta square, red upside-down triangle, and cyan circle refer to a DRG (van D okkum et al.| [2004). to the median of SMGs (Small et al. 20Q% 
[Swinbank et al. 2004), and to MS 1512- cB58 (Te plitz et al. 2000; Baker et al. 2004). respectively. The pink hatched band shows the region occupied b y UV- 
selected BX/BM galaxies at z ~ 2 from lErb et al.i <2006d) . Solid lines represent respectively the fits for the M^-SSFR relations at z ~ 2 (indigo; D addi et al.l 
|2p07a), and those at z — 1 and z — (green and orange, respectively; lElbaz et aL.2007i) . Dotted and dashed lines associated with solid lines represent the la 
scatter of the distributions. 
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Figure 16. Model tracks from PEGASE.2 overplotted on the M*-Z relations presented in Figure fTS] Top and bottom panels refer to closed-box and infall models, 
respectively. Red and light-blue dashed Hnes correspond to = 0.1 and 5 Gyr models, respectively. Three parallel lines for each color represent the initial total 
gas mass as indicated at the top of each panel. The numbers near the small open circles along the tracks indicate galaxy ages in Gyr. The corresponding age for 
each open circle is the same for the lines with the same color. 
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Figure 17. Model tracks from PEGASE.2 overplotted on Figure[T5] Top and bottom panels show tracks from closed-box and infall models, respectively. Red 
and light-blue dashed lines con'espond to r^f = 0. 1 and 5 Gyr, respectively. Three parallel line sets of each color represent the initial total gas mass of as indicated. 
Small open circles and texts along the tracks indicate galaxy ages in Gyr Note that the metallicity with given Tsf and age is independent of Mtotal- 
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1 1 24 34.4 


-21 45 49.4 


2.14 


25.62 


24.77 


24.30 


23.89 


22.89 


21.97 


1.73 


1.92 


10.76 


0.46 


D3a 6397 


11 25 10.5 


-21 45 06.1 


2.00 


24.14 


23.32 


22.78 


22.43 


21.27 


20.48 


1.71 


1.95 


11.29 


0.45 


D3a 7429 


11 25 12.9 


-21 43 30.1 




24.74 


24.14 


23.52 


23.39 


22.27 


21.46 


1.35 


1.93 


10.89 


0.36 


D3a 8608 


1 1 25 06.3 


-21 41 52.7 


1.53 


23.29 


22.99 


22.83 


22.71 


22.09 


21.50 


0.58 


1.21 


10.72 


0.17 


D3a 11391 


1 1 24 54.7 


-21 34 13.3 




23.37 


23.12 


22.73 


22.91 


21.80 


21.06 


0.46 


1.85 


11.16 


0.14 


D3a 12556 


11 25 11.6 


-21 35 48.6 


1.58 


23.44 


22.96 


22.76 


22.59 


21.60 


21.13 


0.85 


1.46 


10.92 


0.24 


Dad 1901 


14 49 41.4 


8 59 50.5 


1.60 


23.56 


23.02 


22.90 


22.48 




20.57 


1.08 


1.91 


11.24 


0.30 


Dad 2426 


14 49 29.0 


8 51 52.9 


2.36 


25.07 


23.92 


23.56 


23.04 




20.75 


2.03 


2.29 


11.25 


0.53 


Dad 2426-b'' 




























Dad 3551 


14 48 59.6 


8 52 06.6 


1.60 


23.54 


22.98 


22.93 


22.74 




20.87 


0.80 


1.87 


11.10 


0.21 


Ha non-detection 


D3a 6004 


11 25 03.8 


-21 45 32.7 


2.39 


25.88 


24.69 


24.31 


23.93 


22.79 


20.96 


1.95 


2.97 


11.32 


0.51 


D3a 6048 


11 25 23.1 


-21 45 32.2 


1.41 


24.33 


23.88 


23.28 


23.01 


21.99 


20.94 


1.32 


2.07 


11.13 


0.36 


D3a7182 


11 25 12.4 


-21 43 48.1 




26.74 






24.68 


23.41 


21.99 


2.06 


2.69 


10.84 


0.54 


D3a 8249 


11 25 15.9 


-21 42 23.3 


2.06 


25.25 


24.38 


23.47 


22.98 


21.54 


20.45 


2.27 


2.53 


11.42 


0.59 


D3a 12153"^ 


1 1 24 40.7 


-21 36 19.2 


2.34 


25.04 


24.58 


24.29 


24.24 




21.82 


0.80 


2.42 


10.85 


0.23 


D3a 13557 


1 1 25 04.9 


-21 37 09.1 


2.40 


26.15 


25.30 


24.53 


24.21 


22.75 


21.58 


1.94 


2.63 


10.99 


0.51 


D3a 13600 


11 25 05.9 


-21 37 13.1 


2.27 


24.71 


24.28 


24.02 


23.89 


22.78 


22.08 


0.82 


1.81 


10.62 


0.23 


D3a 14009"^ 


11 24 34.4 


-21 38 34.1 


2.43 


24.59 


23.65 


23.56 


23.37 




21.81 


1.22 


1.56 


10.67 


0.33 


Dad 759 


14 50 11.4 


9 07 51.3 




26.93 


25.21 


25.16 


24.36 




20.94 


2.57 


3.42 


11.42 


0.67 


Dad 1250 


14 49 59.2 


8 56 18.0 




25.70 


24.88 


24.54 


24.08 




21.01 


1.62 


3.07 


11.32 


0.43 


Dad 2079 


14 49 37.3 


8 58 16.6 


0.73'' 


25.07 


24.12 


23.45 


23.19 




20.95 


1.88 


2.24 


11.16 


0.50 


Dad 2742 


14 49 20.5 


8 50 52.4 


2.12 


24.83 


24.02 


23.67 


22.95 




20.84 


1.88 


2.11 


11.18 


0.50 


Dad 3882 


14 48 49.0 


8 53 15.1 


2.25 


25.12 


24.01 


23.76 


23.23 




20.99 


1.89 


2.24 


11.15 


0.50 


Dad 3977 


14 48 45.8 


8 57 45.8 




25.00 


23.85 


23.50 


22.89 




20.23 


2.11 


2.66 


11.54 


0.55 


Dad 4008 


14 48 44.9 


8 53 46.9 


2.22 


27.31 


25.97 


25.84 


25.03 




21.05 


2.28 


3.98 


11.50 


0.59 



Note. — All magnitudes are shown in the AB system and are measured within a 2" aperture |KQ6|) . 
^ Derived from the S^A'-cahbration of lDaddi et al] |2QQ4ai ). See Equations [T] and (2) in PTTI 
^ This object is not included in our ^T-selected catalog, but it was serendipitously found in the FoV of SINFONI. 
^ This object is not included in our current ST-selected catalog. Quoted values are from an older version of the catalog. 

This object was considered at 2 — 2 based on its photometric redshift, but it was later spectroscopically confirmed to be at c = 0.73. 
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Table 2 

Observational Information 



ID 


Date 


Instrument 


Grism 


Exposure 










(seconds) 


Ha detection 


UJa jzo / 


luuj April Ij 


^IWrUJNl 


H+K 


3600 


D3a 4624 


2005 April 16 


SINFONP 


H+K 


3600 


D3a 4654 


2005 April 16 


SlNFONl 


H+K 


3600 


D3a4751 


2005 April 16 


SINFONI 


H+K 


3600 


D3a5814 


2005 April 14 


SINFONI 


H+K 


5400 


D3a 6397 


2005 April 15 


SINFONI 


H+K 


3600 


D3a7429 


2005 April 15 


SINFONI 


H+K 


3600 


D3a 8608 


2005 April 25 


OHS 


JH 


8400 


D3a 11391 


2005 April 16 


SINFONI^ 


H+K 


1800 


D3a 12556 


2005 April 25 


OHS 


JH 


3600'' 




2005 April 26 


OHS 


JH 


3600 


Dad 1901 


2004 May 1 


OHS 


JH 


8000 


Dad 2426 


2004 May 6 


CISCO 


wK 


6000 




2005 April 14 


SINFONI 


H+K 


7200 


Dad 3551 


2005 April 24 


OHS 


JH 


14400 




2005 April 26 


OHS 


JH 


7200"^ 


Ha non-detection 


D3a 6004 


2005 April 16 


SINFONI 


H+K 


3600 


D3a 6048 


2005 April 24 


OHS 


JH 


10800 


D3a7182 


2005 April 16 


SINFONI 


H + K'' 


3600 


D3a 8249 


2005 April 16 


SINFONI 


H+K 


3600 


D3a 12153 


2005 April 14 


SINFONI 


H+K 


3600 


D3a 13557'i 


2004 May 6 


CISCO 


wK 


7000 




2005 April 30 


CISCO 


wK 


7000 


D3a 13600'' 


2004 May 6 


CISCO 


wK 


7000 


D3a 14009 


2005 April 30 


CISCO 


wK 


5000 


Dad 759 


2005 April 16 


SINFONI 


H+K^ 


1800 


Dad 1250 


2005 April 16 


SINFONI 


H+K^ 


1800 


Dad 2079 


2004 May 6 


CISCO 


wK 


3000 


Dad 2742 


2005 April 15 


SINFONI 


H+K 


7200 


Dad 3882 


2005 April 16 


SINFONI 


H+K 


1800 


Dad 3977 


2005 April 25 


OHS 


JH 


2400 


Dad 4008 


2005 April 25 


OHS 


JH 


3600 



^ AO-module was used. 



" This exposure was not used because it does not improve S/N ratio due to 
worse seeing. 

' This exposure was not used because it does not improve S/N ratio due to 
reduced signal by cloud passage. 

^ D3a 13557 and D3a 13600 were supposed to be in the same slit, but the 
objects were not correctly placed. 

The target acquisition was failed and the object was not correctly placed in 
the slit. 



Table 3 

Emission Line Properties of Ha and [N II] A6583 



ID 




/(Ha) 


EW.-estCHa)"-'' 


/([N II]) 


EW,-es,([N ll])" 


[N ii]/Ha 








(10""ergs s^'cm"-) 


(A) 


(10""ergs s^'cm"-) 


(A) 




km s"' 


















Dad 1901 


1.600 


15± 1 


38 ±8 


7.9 ±0.7 


20 ±4 


0.53 ±0.06 




Dad 3551 


1.603 


17± 1 


55 ± 12 


1.6 ±0.9 


5±3 


0.09 ±0.05 




D3a 8608 


1.528 


16±2 


86 ±20 


11 ± 1 


60 ± 13 


0.70 ±0.1 1 






1.588 


5U±3 


184 ± 38 


1 1 ± 3 


4U ± 14 


O.ii ± O.Oo 




Cisco 


















2.397 


9.0ib 1.8 


22 ± 6 


3.4 ± 2.7 


9 ± 7 


n :iQ _|_ n -1 1 
\j . J o zl \} ^ 




SiNFONI 
















Dad 2426-b 


1.772 


31±1 


> 194 


12 ± 1 


> 76 


0.39 ±0.04 


245 


D3a 3287 


2.205 


9.2 ±2.0 


89 ±26 


< 2.0 


< 19 


<0.22 


118 


D3a 4626 


1.586 


4.4 ±0.3 


22±6 


< 3.0 


< 15 


<0.68 


212 


D3a 4654 


1.551 


33±2 


52±11 


21±1 


33±7 


0.63 ±0.05 


170 


D3a4751 


2.266 


21±2 


165 ±37 


4.1±1.6 


33± 14 


0.20 ±0.08 


74:'! 


D3a5814 


2.141 


39±3 


316±68 


15±3 


120 ±34 


0.38 ±0.08 


169 


D3a 6397 


1.513 


47±5 


99 ±23 


16±5 


34± 13 


0.34±0.11 


265 


D3a 7429 


1.694 


17±4 


96 ±29 


6.0 ±3.5 


32 ±20 


0.34±0.21 


125 


D3a 11391 narrow 


1.774 


13±1 


46 ±10 








187 


broad 




35 ±2 


123 ±26 








1040 



CD 

n 

o 
n 

Note. — [N ii] stands for [N ii] A6583. The object IDs with slanted fonts are considered to be AGN-dominated .vfiz^^s judged from a large [N ii]/Ha ratio ^ 
(D3a-8608) or a broad-hne component {D3a-1 1391). t§ 
^ Equivalent widths are calculated by assuming the continuum flux from the best-fit SED. For Dad-2426-b, K^i = 21.5 |K06|) is used for upper limit of the ^ 
continuum ^ 
^ Ha equivalent widths are con'ected for underlying stellar Ha absorption calculated from the best-fit SED. For Dad-2426-b, the average of the Ha ^ 
absorption line equivalent widths of the other objects, 4.4 A, is assumed. ^ 
The velocity dispersions are con'ected for the instrumental resolution, tJinsirnmcni = 85 km s~^ ^ 
This is smaller than instrumental resolution, hence the uncertainty is large. ^ 

? 

to 
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Table 4 

SED Fitting Results 





(mag) 


10g(,M*SED / ™ .) 


Ohs 






Dad 1Q01 


0.25 


11.23 


Dad 'iS'il 


0.20 


11.08 


n'^a 8608 


0.25 


10.20 




0.20 


10.73 


Cisco 






Dad 2426 


0.65 


11.63 


SiNFONI 






Dad 2426-b 






D3a 3287 


0.15 


10.75 


D3a 4626 


0.52 


10.94 


D3a 4654 


0.45 


11.43 


D3a 4751 


0.25 


10.72 


D3a5814 


0.60 


10.90 


D3a 6397 


0.55 


11.14 


D3a 7429 


0.40 


10.72 


D3a 11391 


0.10 


11.10 



^ E(B-V) from the stellar continuum. 
Derived by SED fitting. 



Table 5 

Extinctions and Hq Star-formation Rates 



ID 




4 gas 


Ana 




SFRraw^ 


SFRcorr*^ 




(mag) 


(mag) 


(mag) 


(10'*^ ergs s"') 


(Mq yr-1) 


(M0 yr-l) 


Ohs 














Dad 1901 


1.0 


2.1 


1.7 


2.2 ±0.2 


20±1 


98±7 


Dad 3551 


0.6 


1.4 


1.1 


2.7 ±0.2 


24±1 


87±5 


D3a 8608 


1.0 


2.1 


1.7 


2.2 ±0.2 


19±2 


92 ±12 


D3a 12556 


0.8 


1.7 


1.4 


8.1 ±0.6 


66±4 


246 ±15 


Cisco 














Dad 2426 


2.6 


5.1 


4.2 


3.3 ±0.8 


33±7 


1540 ±320 


SiNFONI 














Dad 2426-b 








6.1 ±0.1 


53±2 




D3a3287 


0.6 


1.4 


1.1 


3.2 ±0.7 


28±6 


77 ±17 


D3a 4626 


2.1 


4.1 


3.4 


0.6 ±0.3 


5.7 ±0.4 


129 ±8 


D3a 4654 


1.8 


3.6 


3.0 


4.6 ±0.3 


41±3 


622 ±39 


D3a 4751 


1.0 


2.1 


1.7 


7.9 ±0.6 


66±6 


323 ±31 


D3a5814 


2.4 


4.7 


3.9 


13±1.0 


110±8 


3880 ± 290 


D3a 6397 


2.2 


4.3 


3.6 


6.5 ±0.7 


55 ±6 


1460 ±160 


D3a 7429 


1.6 


3.2 


2.6 


3.3±0.7 


28 ±6 


315±71 


D3a 1139P 


0.4 


1.0 


0.8 


2.5 ±0.3 


23±2 


48±4 



" Values not corrected for extinction. 
Extinction corrected values. 

Only the contribution from the narrow-line component is listed and the narrow-line is assumed to 
come from star-forming regions. 
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Table 6 

Gas-phase oxygen 
abundances 



ID 12+log(0/H) 



Ohs 
Dad 1901 
Dad 3551 
D3a 8608 
D3a 12556 

Cisco 
Dad 2426 

SiNFONI 
Dad 2426-b 
D3a 3287 
D3a 4626 
D3a 4654 
D3a4751 
D3a 5814 
D3a 6397 
D3a7429 
D3a 11391 



9.21 ±0.09 
8.55±0.21 
(9.53 ±0.24) 
8.84 ±0.08 



9.02 ±0.34 



9.03 ± 0.04 
< 8.84 
<9.48 
9.38 ±0.10 
8.81 ±0.12 
9.02 ±0.09 
8.97 ±0.12 
8.97 ±0.22 



Note. — D3a 8608 which has 
[N ll]/Ha > 0.7 is considered to 
host an AGN. Thus the derived 
metallicity is listed within paren- 
theses. 



Table 7 

Dynamical masses for 
SINFONI objects 



ID 


Iog(Mdy„/Mo) 


Dad 2426-b 


11.6 


D3a 3287 


11. 1 


D3a4626 


11.5 


D3a 4654 


II. 3 


D3a475I 


10.6 


D3a5814 


11.4 


D3a 6397 


II.7 


D3a 7429 


II.O 


D3a II39I 


1 1.4'' 



''' Calculated from the veloc- 
ity dispersion of the narrow-line 
component. 



